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Abstract 

We show that for a class of dynamical systems, Hamiltonian with respect to 
three distinct Poisson brackets {Pq, Pi, P2), separation coordinates are provided 
by the common roots of a set of bivariate polynomials. These polynomials, which 
generalise those considered by E. Sklyanin in his algebro-geometric approach, 
are obtained from the knowledge of: (i) a common Casimir function for the two 
Poisson pencils (Pi — XPq) and (P2 — fJ-Po)', (ii) a suitable set of vector fields, 
preserving Pq but transversal to its symplectic leaves. The frameworks is applied 
to Lax equations with spectral parameter, for which not only it establishes a 
theoretical link between the separation techniques of Sklyanin and of Magri, but 
also provides a more efficient "inverse" procedure to obtain separation variables, 
not involving the extraction of roots. 

1 Introduction 

The relationship between the Liouville integrability of a Hamiltonian system and the 
existence of a second conserved Poisson bracket (or "hamiltonian structure") in its phase 
space, first discovered by Magri |Q, has been thoroughly investigated in the past years. 
Bihamiltonian structures underlying all classical examples of integrable systems (both 
finite and infinite-dimensional) have been described by several authors, and almost all 
the relevant properties connected to integrability have been reinterpreted in terms of the 
geometry of bihamiltonian manifolds and vector fields. Recently, the classical problem of 
characterizing separable hamiltonians (i.e. those for which the Hamilton-Jacobi equation 
can be solved by separation of variables in a suitable system of canonical coordinates) 
has been also translated in the language of bihamiltonian geometry ^ . 

A question which has not yet received a complete answer concerns the link between 
the bihamiltonian framework and the algebro-geometric methods of solution based on 
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the isospectrality property of Lax equations [Q. Although it is possible to intro- 
duce bihamiltonian structures which naturally lead to Lax equations with a spectral 
parameter |]TT| |2T|, the role of the characteristic equation for the Lax operator (the 



" spectral curve" of the algebro-geometric approach) has not been clarified so far in the 
bihamiltonian perspective. 

The present work adds new elements in view of a connection between multihamilto- 
nian structures, existence of separation coordinates and spectral curves, starting from 
an apparently marginal observation: some well-known integrable systems allow two 
distinct bihamiltonian descriptions, independently described by different authors and 
apparently unrelated (in spite of having one Poisson bracket in common). In this in- 
troductory section, we shall recall some relevant facts using the simplest example of 
such "trihamiltonian systems", namely the generalized Euler-Poinsot rigid body. To 
motivate the reader to follow us through an exercise which could seem of little practical 
interest, let us anticipate that the interplay of the three Poisson structures leads to a 
new role played by the characteristic determinant of the Lax matrix, and this fact may 
eventually clarify the connection between Sklyanin's algebro-geometrical construction 



of separation variables |12| and the bihamiltonian method recently proposed in 00. 

Indeed, the occurrence of more than two Poisson brackets on the same manifold is 
not new nor surprising by itself, and in some cases it is even a structural property, as for 



the so-called "Lie-Poisson pencils" described in |TT|; in the sequel, we discuss the differ- 
ence between such known cases of multihamiltonian structures and the trihamiltonian 
structure that we are presently considering. 

The simplest (nontrivial) example of Lax equation with spectral parameter is pro- 
vided by the dynamics of a rigid body about a fixed point, in the absence of external 
forces (Poinsot rigid body). In the body reference frame, the motion is described by the 
Euler-Poisson-Lax equation 

^^HM."], (1.1) 

where M and VL are the skew-symmetric 3x3 matrices representing the angular momen- 
tum and the angular velocity, respectively, in the body reference frame. A straightfor- 
ward consequence of (|1.1| ) is that the trace of any power of the matrix M is a constant 
of motion: J^Tr(M'^') = 0. Generalizing the system to M G so(r), one obtains in this 
way at most | independent constants of motion (if r is even, or for r odd), which 
for r > 3 would not be enough to meet Liouville's integrability condition. 

Assuming (Ji, l2,h) to be the eigenvalues of the inertia tensor, one can introduce the 
diagonal matrix with diagonal elements ( -^^+^^+13 ^ h-h+h ^ h+h-h y^ ^j-^g linear relation 
between M and Q can then be written in the following form: 

M = m + VLJ] (1.2) 
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Manakov has observed that the Euler-Poisson equation ( [L . 1| ) and the inertia map 



(|1.2|) can be put together into a single Lax equation for a new Lax pair dependent on a 
formal parameter A, 



2^ 



d(M + A J , , , , , 

^ — = [M + \j',n + \J]. (1.3) 

The insertion of the parameter A into the Lax equation leads to a wider number of 
constants of motion, the Manakov integrals. We denote them by fi, according to the 
following convention: 

/f = ij^^A^ + X^/fA'=-\ (1.4) 

1=1 

For M G so(r), the functions fi{M) vanish identically for i odd; the odd Manakov 
functions are however relevant for the "generalized Euler-Poinsot system", with M G 
g[(r), that we shall consider in the sequel. 

As is well known, the equation ( |1 . 1|) is Hamiltonian with respect to the Lie-Poisson 
bracket, defined on so(r) through the ad-invariant scalar product 

{A,B) =Tr{A- B). (1.5) 

More precisely, given any function / : so(r) R, one defines the gradient at a point 
M as the matrix V/ G so(r) such that / = {df,M) = [Mj'Vf); then, for any pair of 
functions, 

{f,g},, = iM,[Vf,Vg]) (1.6) 



is a Poisson bracket |2^. The Lie-Poisson bracket ( |l.(j| ) is degenerate: an ad-invariant 



function f{M) is in involution with any other function, i.e. is a Casimir function for 
the bracket ( |1.6|) . The Casimir functions, which include the traces of the powers of M, 
are automatically constants of motion, but they are irrelevant as far as the Liouville 
integrability of the system (with the Lie-Poisson bracket) is concerned. Therefore, the 
integrability of the system actually relies on the existence of the other Manakov first 
integrals (among which one can find enough independent functions). 

In 1996, Morosi and Pizzocchero |]T3[ introduced a second Poisson bracket on so(r). 



defined as follows: let A G 0l(r) a fixed matrix (for the Euler-Poinsot case, A = J^; 
notice that A needs not to belong to 5o(r)). With the same definition of scalar product 
and gradient as above, one sets 

{/, g}^.^^ = iM,\/f-A-Vg-Vg-A- V/). (1.7) 

One can check that the vector field generated by the Manakov functions through the 
Poisson structure (|1.7| ) are exactly the same as those generated through the Lie-Poisson 
structure ( |1.6| ), up to a rearrangement in the correspondence between hamiltonians 
and vector fields. For instance, the physical hamiltonian generating the Euler-Poinsot 
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dynamics through the Lie-Poisson bracket is hi = ^Tr{QM) = Tr(yi?J), while the 
hamiltonian of the same vector field through the Morosi-Pizzocchero bracket is /i2 = 
— Yl'r{VLJ~^ M J~^) . To simplify the notation, let us denote by Pi and P2 the Poisson 
tensors associated respectively to the brackets ( |1.6|) and (|L 



{f,g},, = {df,Pidg), {f,g},,, = {df,P2dg). (1.8) 

Denoting by Xi the vector field over 50 (r) corresponding to equation ( |1.1| ), the relation 
Pidhi = P2dh2 is depicted by the diagram 



Xi 

hi /i2 



;i.9) 



p p 

where /i — > X is an abbreviation for dh 1 — > X, a convention that we shall use in analo- 
gous diagrams throughout this article. The diagram ( p.. 91 ) is nothing but the elementary 
block of the Lenard-Magri recursion generating a whole family of quadratic first inte- 
grals hi (known as Miscenko functions) , and the corresponding symmetry vector fields 
Xf. 

Xi X2 ■ ■ ■ Xk 

Pl^ >\P2 Pl^ XP2 Pp^ XP2 Pp^ XP2 



hi h2 hs hk hk+i (1-10) 

The Manakov first integrals can be generated by the same recursion procedure. Setting 
A = J2, one has V/f ^ = (M + XA)''-\ then 

Pidft = [M, V/f] = AVft'M - MVft'A = P2dft' (1.11) 
which correspond to Lenard-Magri diagrams starting with the Pi-Casimir functions 



^k ^k 
py Pyr xPi Py^ Pyf 

fk rk+i 
Jk Jk Jk 

Notice that all the functions iteratively generated by Lenard-Magri recursion rela- 
tions are automatically in involution with respect to both Poisson tensors Pi and P2- 
The (elementary) proof of this fact will be recalled in the next section. Thanks to 
this property of bihamiltonian vector fields, one does not need to prove separately the 
involutivity of the first integrals of Manakov, and the complete integrability of the gen- 
eralized Euler-Poinsot system is simply assessed by computing how many independent 
first integrals can be found in this way. 



4 



All the statements above hold valid if one extends the equation ( |1 . 1| ) to M G Qi{r). 
Both the Lie-Poisson bracket and the Morosi-Pizzocchero bracket can be introduced 
in g[(r) using the same definitions ( p,.6D and ( |1 . 7] ) . The Morosi-Pizzocchero bracket 
is defined in terms of the matrix product (not of the commutator) and therefore is 
defined even more naturally on g[(r): it reduces on so(r) by orthogonal projection with 
respect to the scalar product (|1.5| ), provided the matrix A is symmetric. Thus, for 
the Lax matrix L{X) = AX + M the complete family of Manakov constants of motion 
can be found by the recursion procedure, which also ensures their mutual involutivity. 
Whenever A is symmetric (and positive), the dynamical system defined by ( |1 . 1|) in Ql{r) 
is a proper extension of the original Euler-Poinsot system. The fiows of the original 
model are those associated to the even Manakov functions (these fiows are tangent to 
3o(r)), while the other fiows of the enlarged system, generated by the odd Manakov 
functions, are orthogonal to so(r). 

In the larger phase space Ql{r), however, one can obtain the full set of first integrals 
by another Lenard-Magri recursion, relative to a different bihamiltonian pair. The new 
Poisson bracket depends, as for (|1.7| ), on the choice of the matrix A: 

{f,g}, = {A,[Vg,Vf]). (1.12) 

From now on, let us denote by Pq the Poisson tensor associated with this bracket. The 
Manakov functions are in bihamiltonian recursion also with respect to the pair [Pq, Pi), 
but the sequences are arranged in a different way: 

Pid/f = [M, V/f ] = A] = Porf/ti- (1.13) 

Thus, each integer power of L{X) corresponds to a single finite Lenard-Magri sequence, 
starting from a Casimir function for Pq and ending with a Casimir function for Pi: 





ly^ ly ry 

fk fk fk 

Jk J2 Jl 

(1.14) 

A disadvantage of the new bihamiltonian structure (Pq, Pi) is that it cannot be reduced 
(by restriction or by orthogonal projection) to so(r). On the other hand, (Pq, Pi) leads 
naturally to the Lax equation with spectral parameter ( |1.3| ), which on the contrary is 
rather difiicult to derive from the former pair (Pi, P2). To show this, we need to reexpress 
the Lenard-Magri recursion relations ( |1.14| ) in the language of Poisson pencils. 

Given a pair of Poisson tensors {P,Q) on a manifold A4, assume that the A- 
dependent bracket 

{/, 9}. -.a = {/> 9}p - Hf, 9}o = (df, (P - \Q)dg); (1.15) 
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be a Poisson bracket, i.e. fulfill the Jacobi identity for any A; in this case, P and Q are 
said to be compatible; {Ai, P,Q) becomes a hihamiltonian manifold (or PQ-manifold, 
following [Q), and one refers to the bracket ( |1.15D as to its Poisson pencil. It is immediate 
to see that, given a sequence of functions {/i}j=o,...,Ar such that 

= Pdfo 

Qdfo = Pdh 

\ (1.16) 
Qdfk = Pdfk+i 

QdfN = 0, 

then the polynomial in A defined by 

/a = /o + /iA + --- + /,vA^ (1.17) 

is a Casimir function of the Poisson pencil, i.e. for any A 

{P-\Q)df^ = Q (1.18) 

(the differential of fx is taken with respect to the coordinates on A^, A being regarded 
as a parameter). Conversely, given a A-polynomial function fulfilling ( |1.18|) , its coeffi- 



cients obey the Magri-Lenard recursion according to ( |1.16| ) and generate a sequence of 
commuting bihamiltonian vector fields. 

In the next section, we will recall the proof of the following relevant property, that 
we shall extensively use. Let gx be a second Casimir function of the same Poisson 
pencil: then, not only its coefficients Qk are in involution among themselves, but they 
also Poisson-commute with all the coefficients fk of the other Casimir function fx- 

Given a polynomial Casimir function /a, each bihamiltonian vector field of the as- 
sociated Lenard-Magri hierarchy Xk = Pdfk = Qdfk-i can be also represented by a 
Hamilton equation with spectral parameter. Having set fx = YliLofi^^^ consider for 
each positive integer k < N the polynomial 

/f ^ foX' + f,X''' + ■■■ + fk-iX + fk; (1.19) 



taking into account (|1.16|) it is easy to see that 

Xk = {P - XQ) dfi'K (1.20) 

This formula holds true for formal power series (A^ = oo); if the Casimir function fx is 
instead expanded in Laurent series, fx = Xli^o fi'^~^^ then the polynomial f^'' is easily 
obtained by multiplication by A'^ and truncation to the nonnegative powers: 

/f = (AVa)^ . (1.21) 
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We are now ready to derive the Manakov equation ( |1.3D as a Hamilton equation with 
spectral parameter for the Poisson pencil (Pi — APq) on g[(r). 

In fact, it is easy to see that the trace of any power of the Lax matrix AX + M is a 



Casimir function of the Poisson pencil (Pi — APq): by definition ( |1.6| , p. .121 ) 



(Pi - APo) df = [AX + M, V/]; (1.22) 

as already seen, for ^ = iTr(AA+M)^ one has Vf[^^ = (AA + M)'=-\ which obviously 
commutes with AX + M. The same happens for the Laurent series expansion of the 
trace of any half-integer power of A + MX~^. On account of (|1.20| ), all the vector fields 
generated by the coefficients of these Casimir functions (which all mutually commute, 
by the property mentioned above) correspond to Lax equations with spectral parameter: 

X, = (Pi - APo) rf/f = [AX + M, V/f ]. (L23) 

In particular, the Manakov equation ( |1.3| ) corresponds to the first vector field of the 
hierarchy associated to the Casimir function fx = |Tr(yl + MA~^)^/^, for which V/}^"* = 
A^/^A + f2. No comparably simple and natural connection exists between the other 
Poisson pencil (P2 — APi) and the Lax-Manakov form of the equations. This setting can 
be generalized to cover the cases of more general Lax matrices with spectral parameter 
on g[(r), of the form L{X) = AX^ + MiA'^~^ + ■ ■ ■ + M^. The general framework is 



described in |]TT|, and will be partially recalled in section (4) below. 

So far, we have simply reviewed some results already present in the literature. Now, 
some questions arise naturally. We have described a dynamical system which is bihamil- 
tonian with respect to two independent PQ structures, (Pq, Pi) and (Pi,P2); is that a 
pure accident, or it is a common situation? 

One can check by explicit computation that the Poisson tensors Pq and P2 are not 
only separately compatible with Pi, but also compatible with each other (that is not ob- 
vious, as compatibility is not a transitive relation). Does it make any sense to introduce 
the notion of a trihamiltonian structure (Po,Pi,P2)? Would it carry any additional 
information not already contained in either one of the PQ structures, each of which 
already allows to characterize completely the dynamical system and its symmetries? 



The vector fields ( |1.23|) on g[(r) are indeed trihamiltonian. The full set of hamilto- 



nians and vector fields generated by the traces of integers powers of AX + M, fit into a 
single "planar" diagram (as was first pointed out by M. Ugagha [0]), which could be 
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regarded as the "trihamiltonian version" of the Lenard-Magri "hnear" diagrams ( |1.14| ) : 





il 

h 

[M, A\ 

X^i Pyf x^i -fy 

fl fi 

sPi XPi fy^ XPi ryf 

fx fi fi 



P2 

(1-24) 

We have seen above that, for a PQ structure, any hnear recursion starting from 
a Casimir function of Q and ending with a Casimir function of P corresponds to the 
existence of a A-polynomial Casimir function of the Poisson pencil (P — XQ) . Can one 
find a "generating polynomial" for the full trihamiltonian recursion? The answer is yes: 
as we shall see in detail in section (2), if one considers two compatible Poisson pencils 
{Pi — XPq) and (P2— A^-Po); one can define a common Casimir function of the two pencils 
to be a bivariate polynomial /a^ = Yl ^] ^''f^^ such that 

(Pi - APo) c^/am = 

(P2 - /iPo) rf/A^ = (1.25) 

for any value of {X, fi): then its coefficients /;,*• fulfill the recursion relations represented 
in the diagram (|1.24|) . Later on we will explain why the construction of two Poisson 
pencils, each one with its own spectral parameter, is here more fruitful than introducing 
a two-parameter pencil like (Pq — APi — /iP2). 

Up to this point, the reader might still regard the idea of trihamiltonian structures as 
an artifact of purely academic interest, a mere "variation on the theme" of bihamiltonian 
structures. Two results, presented in this article, suggest that the subject is worth 
investigating further. 

First, the trihamiltonian structure associate to equation ( |1.3|) can be generalized, 
in quite a nontrivial way, to Lax equations for matrices of the form L{X) = AX"^ + 
MiA""'^ + ■ ■ • + M„, which include several interesting systems such as the Lagrange top 
[P!7| and the finite-dimensional Dubrovin-Novikov reductions of the Gel'fand-Dickey 
soliton hierarchies |jl8|. Indeed, the generalization of the pencil (Pi — APq) to the direct 
sum of n copies of Ql{r) was already described in [|n[], but to our knowledge it is an 
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entirely new result that also the Morosi-Pizzocchero bracket is a particular case of a 
more general structure existing on gt(r)'^, a Poisson tensor P2 which turns out to be 
quadratically dependent on the dynamical variables Mj, apart for the linear case n = 1 
already discussed. 

The second striking fact is that for these trihamiltonian structures on 0[(r)" there 
always exists a common Casimir function of the two pencils, which (for a generic choice of 
the matrix A) has the property that its coefficients form a maximal set of independent 
hamiltonians in involution (we stress that, in contrast, the recursion diagram for the 
traces of the powers of the Lax matrix includes infinitely many hamiltonians, and one 
has to single out a finite subset of independent first integrals). This miraculous Casimir 
function is nothing but the characteristic determinant of the Lax matrix, 

/A^ = det|L(A)-/iI|. (L26) 

The corresponding recursion diagram features a sort of "fundamental molecule" , a "fin- 
gerprint" associated to the trihamiltonian structure (P, Qi,Q2)- For instance, the fol- 
lowing diagram corresponds to the trihamiltonian recursion on the algebra 0t(3): 





1 ry^ 




K 




h 


h 


Xi 


X2 







00 




hi 


\" 




















(L27) 

while ( p..28| ) is the "molecule" of the trihamiltonian structure on 0t(2)^: 



h\ 



Xi 



P2 







h\ 

h 
X2 



hi 



Pi 







(L28) 
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The general form of the "fundamental molecule" for g[(r)" is given in section (4) as 
fig.l. 

Although it is a well known fact that the coefficients of the characteristic polynomial 
are in involution with respect to the usual Lie-Poisson bracket, in the bihamiltonian 
framework there was no apparent reason to introduce a bivariate polynomial fx^ in 
connection with the Lenard-Magri recursion. For a trihamiltonian structure, instead, 
it is quite natural to consider this object, and the characteristic polynomial of a Lax 
matrix becomes just a particular case of it, in exactly the same way as Lax equations 
with spectral parameter are a particular case of Hamilton equations, for the appropriate 
Poisson pencil ( p..2(J| ). 



This opens a very interesting perspective. The characteristic equation 

det|L(A) -/ill = 0, (L29) 

regarded as a polynomial equation for (A, fi) G defines the well-known spectral curve, 
i.e. the starting point for the algebro-geometric methods of linearisation [P] [§]. In 
the trihamiltonian framework, as we have seen, the characteristic determinant naturally 
occurs as the fundamental Casimir function of two pencils: yet this does not explain why 
the roots (Aj, /Xj) of eq. ( |1.29| ) should play any role at all. Now comes a third surprise: 
a fairly general construction presented in section (3) shows that the equation fx^ = is 
the keystone for the construction of canonical separation coordinates for trihamiltonian 
systems. 

This result essentially derives from an observation by E. Sklyanin |]T^ . On algebro- 
geometric grounds, Sklyanin has found a "magic recipe" ( "Take the poles of the properly 
normalized Baker- Akhiezer function and the corresponding eigenvalues of the Lax oper- 
ator"), which essentially amounts to finding the common roots of ( |L29|) and of suitable 
minors (or linear combination of minors) of the characteristic matrix L{X) — /xl. In 
the examples considered by Sklyanin, the new variables [Xi, fii) defined in this way turn 
out to be canonical with respect to a suitable Poisson bracket; by direct consequence of 
eq. ( p..29| ), all the hamiltonians defined as the (nonconstant) coefficients of fx^ are then 
separable in the coordinates {Xi,^i). However, Sklyanin himself remarks that "generally 
speaking, there is no guarantee that one obtains the canonical Poisson brackets [..] The 
key words in the above recipe are 'the properly normalized'. The choice of the proper nor- 
malization can be quite nontrivial, and for some integrable models the problem remains 
unsolved". 

Independently of Sklyanin's approach, Magri and his collaborators ||2|||^ have re- 
cently shown that given (i) a PQ structure, (ii) a complete family of commuting hamil- 
tonians defined by the Casimir functions of the Poisson pencil, and (iii) a set of vector 
fields, suitably normalized on the hamiltonians previously introduced, which preserve 
the Poisson tensor P but do not belong to its image (geometrically speaking, they should 
be transversal to the symplectic leaves of P), then one can define by projection (under 
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some additional conditions on the vector fields) a reduced, kernel-free bihamiltonian 
structure; for this new PQ structure, a set of Darboux-Nijenhuis canonical coordinates 
can be obtained by a constructive procedure, and the original hamiltonians (properly 
reduced) turn out to be all simultaneously separable in these coordinates. 

The theoretical interest of both constructions is largely beyond the concrete appli- 
cability of these procedures. As a matter of fact, while Sklyanin's recipe lacks a general, 
theoretically-grounded rule to find the key element (the normalization of the BA func- 
tion, or equivalently the proper linear combination of minors of the Lax matrix which 
should vanish), in Magri's theory there is no practical recipe to construct systematically 
sets of transversal vector fields fulfilling the necessary requirements. In both approaches, 
moreover, the final construction of separation coordinates involves finding the roots of 
polynomial equations, which even for rather simple examples turn out to be of order 
higher than three. 

As we show in this article, Magri's procedure can be adapted to the trihamiltonian 
setup, without loosing its geometric elegance, and actually making the theory even sim- 
pler and more symmetric (although less general). In this framework, the central role of 
the "generalized spectral equation" /a^ = becomes clear. Moreover, for the trihamilto- 
nian structures that we are introducing on the spaces 0t(r)", we have found a systematic 
way to produce explicitly the required transversal vector fields, and we will show that 
the associated Darboux-Nijenhuis coordinates are exactly the roots of suitable combi- 
nations of minors of the characteristic matrix, much alike Sklyanin's coordinates; in this 
way, we provide for this class of systems the missing element in both Sklyanin's and 
Magri's prescriptions for the construction of separation variables. In addition, we show 
that our framework makes available a different strategy, which yields the inverse trans- 
formation (i.e. the matrix elements of the Lax operator as functions of the separation 
variables) by solving only a system of linear algebraic equations, thus bypassing the 
problem of finding roots of higher-order polynomials. 

Let us quote another important remark by Sklyanin |T^: "Separation of variables, 
understood generally enough, could be the most universal tool to solve integrable mod- 
els [...] the standard construction of the action-angle variables from the poles of the 
Baker- Akhiezer function can be interpreted as a variant of separation of variables, and 
moreover, for many particular models it has a direct quantum counterpart" . Therefore, 
a satisfactory hamiltonian setup for Sklyanin's construction is likely to provide a link 
between hamiltonian and algebro-geometric integrability. In this sense, the equation 
fxfi = should deserve some additional interest, as it points towards a generalisation of 
the notion of spectral curve not relying on Lax representations. 

The article is organized as follows: in section (2) we recall, as synthetically as possi- 
ble, some facts about bihamiltonian structures which are necessary for the subsequent 
discussion; then, we present a theoretical setting of our class of trihamiltonian structures. 
In section (3) we discuss the general method of construction of separation variables, 
i.e. the trihamiltonian version of Magri's construction. We present in detail the proofs 
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of some relevant propositions providing the theoretical background for all applications of 
our framework; furthermore, we show how the components of all relevant objects (Pois- 
son structures, common Casimir function, transversal vector fields, etc.) look like in 
Darboux-Nijenhuis coordinates; this will be used in section (4) to reconstruct the coor- 
dinate transformation. The fourth and last section is devoted to the application to Lax 
equations with spectral parameter on 0l(r)"^; here we simply list the "ingredients of the 
recipe" without a general proof; this section is intended to present a concrete outcome 
in just enough detail to motivate the reader to deal with the theoretical construction of 
section (3). 

Throughout the article no attempt is made to present a geometric characterisation, 
or classification, of the trihamiltonian structures possessing the specific features con- 
sidered. In particular, we deliberately avoid to encompass these features into a single 
definition of "special trihamiltonian structure". In our discussion, the basic structure 
involved is sometimes presented as a triple of compatible Poisson tensors, but more of- 
ten as a pair of Possion pencils; the assumptions on this basic structure vary according 
to the context. In sect. (2) we just reconsider the notion of "trihamiltonian recursion" 
associated with a common Casimir function of two Poisson pencils, without imposing 
particular conditions on the Poisson tensors besides their mutual compatibility; in such 
generality, indeed, nothing ensures that a common Casimir function exist at all (we 
will present a counterexample). In sect. (3.1) we show how to define a pair of Nijenhuis 
tensors using a set of vectorfields, which should fulfill a number of requirements: these 
Nijenhuis tensors are related with a particular system of Darboux-Nijenhuis coordi- 
nates, and do not act as recursion operators for the trihamiltonian iteration described 
previously (throughout this subsection, the existence of a common Casimir function is 
irrelevant). In section (3.2) the two objects, i.e. the common Casimir function (that we 
now require to be complete in a suitable sense) and the set of transversal vectorfields, are 
eventually combined together to construct a set of bivariate polynomials Sa{X, //): then, 
a number of additional hypotheses are introduced to obtain the main result, i.e. that the 
common roots of these polynomials are Darboux-Nijenhuis coordinates, in which the 
Hamiltonians occurring as coefficients in the common Casimir function are separated 
(in the sense of Sklyanin). Hence, the set of conditions to be imposed on the basic 
structure depends on whether one looks for the mere existence of iterated trihamilto- 
nian vectorfields, or for a possible explicit construction of separation coordinates. In 
conclusion, we feel that at the present stage of understanding of the matter presented 
herein, one could hardly find a simple and general definition which might be regarded 
as truly fundamental. The aim of the article is rather to display a number of nontrivial 
and rigorous (if not yet complete) arguments in favour of further investigations in this 
direction. 
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2 From bi— to trihamiltonian structures 



2.1 Poisson pencils 

As was already done ( PT^ ) in the introductory section, we shall represent a Poisson 
bracket { ■ , ■ } on a manifold Ai by means of a contravariant antisymmetric tensorfield 
P, according to 

{dg,Pdf) = {f,g}. (2.1) 

The names Poisson structure or hamiltonian structure are equivalently used, as is com- 
monly done, to denote both the tensor P and the algebra of differentiable functions on 
M. with the bilinear operation defined by the corresponding Poisson bracket. Of course, 
a contravariant antisymmetric tensorfield P defines a hamiltonian structure only if the 
bracket ( |2.1D obeys the Jacobi identity; this condition corresponds to a differential iden- 
tity on the components of P. 

In most of our applications, the tensor P will not be of maximal rank; thus, the 
subalgebra of functions which are in involution with any other function may include 
non-constant functions, the Casimir functions. The Casimir functions are constant 
of motion for any hamiltonian vectorfield, i.e. for any vectorfield being the image of a 
closed one-form through the Poisson tensor P. Therefore, any trajectory of any possible 
hamiltonian system on that phase space lies entirely on a common level set of all the 
Casimir functions. Generically, such a level set is a submanifold, the dimension of which 
equals the rank of the Poisson tensor. Upon reduction to any of these submanifolds, the 
Poisson tensor becomes invertible and therefore defines a symplectic structure. For this 
reason, the common level sets (for regular values) of the Casimir functions are called 
symplectic leaves. In contrast with the case of symplectic manifolds, the Lie derivative 
the Poisson tensor can vanish along the flow of a given vectorfield X, 

CxiP) = 0, (2.2) 

without X being even locally hamiltonian: the condition ( |2.2| ) can be fulfilled also by 
vectorfields which do not belong to the image of P, and are therefore transversal to the 
symplectic leaves. Such vectorfields will be called weakly hamiltonian. They play an 
important role in the sequel. 

Bihamiltonian structures were first introduced by Magri in |jl|. 

Defininition 2.1: Two Poisson tensors P and Q on a manifold Ai are 
told to be compatible if any linear combination of the two tensors is again a 
Poisson tensor. 

In such situation one can find vectorfields which are hamiltonian with respect to both 
structures, i.e. bihamiltonian vectorfields. In this article, we borrow from the bihamil- 
tonian theory the following facts: 
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(i) Two hamiltonians are associated to a single bihamiltonian vectors field X = 
Pdh = Qdk. Then, one can define two other vectorfields, namely Qdh and Pdk. In 
some cases, these turn out to be bihamiltonian as well, and the procedure can be iterated 
yielding a Magri-Lenard hierarchy of bihamiltonian vectorfields, as in ( |1.1CI| ). 



(ii) Once a Magri-Lenard hierarchy has been constructed, all the vectorfields be- 
longing to it are mutually commuting, and all their hamiltonians are in involution with 
respect to both P and Q. 

(iii) There are basically two ways to produce such hierarchies: if at least one of the 
Poisson tensors (say, P) is nondegenerate, then one can introduce the recursion operator 
(or Nijenhuis tensor) 

N = Q-P-\ (2.3) 

One can prove 0] that for any bihamiltonian vectors field X, the vectors field NX is 
also bihamiltonian, so the hierarchy can be produced by iterated application of the (1,1) 
tensors field N . Alternatively (for instance, if both Poisson tensors are degenerate), one 
can look for Casimir functions of the Poisson pencil {Q — AP), as already described in 
the Introduction. 

The classical proof of the involutivity property, which is the most relevant to our 
purposes, is so simple and elegant that we reproduce it here (further details can be found 
in 0, ilOi): 



Proposition 2.1: Let fx and g\ be two Casimir functions of the Poisson pencil (Q—XP) . 
Assume that f\ and g\ are expanded in power series in the parameter X of the pencil, 

/a = Ei=o fi^' gx = Ei=o 9i>^'- T^^^ {/i' fk} = {9j, 9k} = {fj, gk] = for all j, k: 
this holds for both brackets { , }p and { , }q associate to P and to Q respectively. 

Proof: The conditions {Q — XP) dfx = and {Q — XP) dgx = are equivalent to 
Pdfi = Qdfi+i and Pdgi = Qdgi+i. Moreover, one should have Qdfo = and Qdg^ = 0, 
i.e. the lowest -order coefficients of both expansions should be Casimir function for Q. 
One can assume j < k, without loss of generality. From the definition ( |2.1| ) {fj, fk}p = 
{fj+i,fk}Q = {fj+i, fk^i}p- Whenever k - j is even, applying repeatedly the equality 
one finds that {fj,fk}p = {fr,fr}p = 0, with r = {k — j)/2; otherwise, after {k — j) 
steps one finds {fj,fk}p = {fk, fj}p, which proves the statement by the antisymmetry 
of the Poisson bracket. The same holds for {gj,gk}p, and for the other bracket { , }g. 
Furthermore, applying the same iterative argument one finds {fj,gk}Q = {fj+k, go}Q, 
and the latter bracket vanishes because go is a Casimir function for Q. This proves that 
{fj^9k}Q = for all j, fc. Since {fj,gk}p = {fj+ugk}Q, one has {fj,gk}p = as well. 
□ 
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2.2 Casimir functions and trihamiltonian vector fields 



Assume that a manifold Ai is endowed with three Poisson tensors Po, Pi and P2, pairwise 
compatible. A natural question is whether a set of vectorfields that are hamiltonian with 
respect to all three structures can be generated by the coefficients of some "generating 
function", analogous to the Casimir function f\ above, and whether the corresponding 
hamiltonians would then be automatically in involution. 

One might believe that the obvious generalization of the setting just described would 
consist in introducing a two-parameter Poisson pencil 

Po - APi - flP2 

and seeking for its Casimir functions. Unfortunately, the coefficients of the Taylor series 
in the two parameters (A, /i) do not fit into any useful recursion relation: from the 
Casimir equation 

00 

(Po - APi - /XP2) d/A^ = for /a,^ = 5^ Iviy/jP 

i,j=0 

one gets the relations 

Pod/i° = 

Podh%^ = PidhP^ i = 0, 1,... 

P^dhl^^ = P^dhi J = 0,1,... 

Podhi^l = Pidhi^^ + P2dh{+^ i,j = 0,1,... 

which neither provide trihamiltonian vector fields nor force the functions // to be in 
involution. 

Let us consider instead a function fx^ that is simultaneously a Casimir function of 
the two distinct pencils (Pi — APq) and (P2 — //Pq): 

(Pi - APo) d/A^ = , (P2 - A^Po) dfx^ = . 

In this case we actually obtain the following relations: 

Pidhi = J = 0,1,... 

Podh{ = Pidhi^, i,j = 0,l,... 

P^dhl = i = 0, 1,... 

Pod/i^ = Psc^/ip i,j = 0,l,... 
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(2.4) 

the vectorfields Pq dh'l = Pi dhi_^i = P2 dhj^^ are clearly trihamiltonian. 

Notice that it is possible to find a common Casimir function which can be (formally) 
expanded in a Taylor series with respect to the two parameters X e fi only if both Pi 
and P2 are degenerated Poisson tensors: in fact, for a fixed power of A, the lowest-order 
coefficient in /i must be a Casimir function of P2, while the lowest-order coefficient in 
A for any fixed power of /i must be a Casimir function of Pi. If, moreover, also Pq is 
degenerate, then it is possible to find Casimir functions which are polynomials in A and 
/i. For such functions the recursion diagram is finite. 

In analogy with the bihamiltonian case, one has: 

Proposition 2.2: Given a common Casimir function = ^ hjX'^fi^ of two compatible 
Poisson pencils Pi — APq and P2 — /iPo, all the coefficients h^ are in mutual involution 
with respect to all three Poisson brackets. 

Proof: For any i,j the functions h^^^ = ^h\\^ and h^^^ = Yl^i^'^^ i-^- the total co- 
efficients of fi^ and fi^ in the expansion of fx^, are Casimir functions of the Poisson 
pencil (Pi — APq); then, by proposition ( p.l| ) all their coefficients hi are in involution 
with respect to both Pq and Pi. On the other hand, the functions h^fl^ = ^ hffi'^ and 
= h'^f^'' are Casimir functions of the other pencil P2 — /iPo, hence h{ are in invo- 
lution also with respect to P2. □ 

We remark that there are other possible ways to extend the bihamiltonian framework 
to the case in which there are more than two compatible Poisson structures. The 
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idea of a trihamiltonian vector field was already considered, for example, in |T9[ and 



161), where the three structures P, Q and 5* were however assumed to produce the 
iteration Pdhi = Qd/ij+i = Sdhi^i. In this approach, the third structure only supplies 
an additional relation which links vectorfields already belonging to the same Magri- 
Lenard hierarchy; in our framework, the third structure acts instead as a bridge linking 
different bihamiltonian hierarchies, and so allows to collect a greater number of function 
in involution in a single objet: the common Casimir function. 

In our setup, the structure Pq seem to play a distinguished role with respect to Pi 
and P2. As a matter of fact, it is easy to figure out how to include in the picture also 
the Poisson pencil built from Pi and P2, but the P1-P2 recursion is already included in 
the diagram (|2.4| ), and introducing a third pencil would be redundant. In the recursion 
diagram, all structures appear on equal footing; on the other hand, in the applications 
that we have in mind there is a distinguished structure, so the "symmetry breaking" 
caused by the choice of two pencils is significant. 

We stress that, given three Poisson structures Pq, Pi and P2, mutually compatible 
and such that the two pencils Pi — APq, P2 — /xPo both admit Casimir functions, a 
common Casimir function as required in Prop. (2.2) may not exist at all. An obvious 
necessary condition is that at each point x of the phase space Ai, and for generic values 
of the spectral parameters (A,/i), the two subspaces ker{Pi — APq) and ker{P2 — /uPq) 
should have a nontrivial intersection in T*Ai. For instance, let us consider the space 
M^, with coordinates {xi, 0:2, 0:3, 0:4, X5}, endowed with the three Poisson tensors: 

0100 
0010 
Pfi = I -1 000 
-1000 
0000 

(00010 
00000 
00001 
-10 00 
0-100 

Although the pencils Pi — APq and P2 — APq admit the Casimir functions X3 + AX4 + X^x^ 
and X2 + fiXi — /i^Xs, respectively, it is easy to check that a common Casimir function 
for both pencils does not exist. 

The simple local geometry of our trihamiltonian structures may be clarified by an 
example. Let us consider the "fundamental molecule" ( 1.27 ). The lowest dimension in 




which this diagram can be realised is 9. In fact, the diagram includes six functions, that 
we assume to be independent. Since the three vectorfields in the diagram commute, 
by Frobenius' theorem there exists a coordinate system in which they coincide with 
coordinate vectorfields: Xi = -£j ior i = 1, 2, 3. The diagram shows that the hamiltonian 
/iq is Po-conjugate to x^, while and hi are Po-conjugate to and x^ respectively. 
The other hamiltonians /i2, h\ and /iq are Casimir functions for Pq. Therefore, the 9 
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functions and h'j should be functionally independent, and locally form a coordinate 
system: let x'^ = Hq, = h\, x^ = x'^ = x^ = h\ and x^ = h"^. It can be read 
directly from the diagram ( |1.27D that in these coordinates the three tensors Pq? -Pi and 
P2 have the following matrix components: 



\ 



1 
1 
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000000 
000000 
000000 
000000 
00000 
00000 
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The fact that three independent Poisson tensor can be simultaneously put in canoni- 
cal form is possible only because they are all degenerate (in the realization of minimal 
dimension, they ought to have the same rank), and their symplectic foliations are dif- 
ferent. Having anticipated that the diagram (|1.27|) corresponds to the trihamiltonian 
structure of 0l(3), we have in fact shown that the latter trihamiltonian space admits 
local multicanonical coordinates. On the other hand, although six of the multicanonical 
coordinates coincide with the coefficients of the characteristic determinant of the Lax 
matrix AA + M, the first three coordinates can be found only upon explicit integration of 
the dynamical system. Analogous considerations hold for the "fundamental molecule" 
of any space g[(n)'^, and actually for any finite trihamiltonian recursion diagram (under 
the assumption that all the hamiltonians are independent, which is generically true for 
Q[{nY). The lowest dimension to accommodate a trihamiltonian structure admitting 
multicanonical coordinates is 4 (the "fundamental molecule" contains just one vectors 
field and three hamiltonians); an example is the algebra 0t(2). The reader can easily find 
out the multicanonical form of the three Poisson tensors generating the 0^(2)^ diagram 

(Pi- 

Another type of coordinates, the separation coordinates, can instead be obtained 
explicitly in an alternative way, which does not require the integration of the vectorfields. 
This is explained in the next section. 



3 Separation of variables 

3.1 Darboux-Nijenhuis coordinates 

In this section we shall adapt to the trihamiltonian framework the construction of sep- 
aration variables proposed by Falqui, Magri and Pedroni in 
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The basic notion involved in their construction is the definition of Darboux-Nijenhuis 
coordinates [0. Consider a bihamiltonian structure PQ on a manifold Ai, with 
dim(A^) = 2m, such that at least one of the Poisson tensors (say, P) is nondegenerate. 
Let be the recursion operator defined by ( p.3| ); following Magri [Q, we say that A4 is 
endowed with a Poisson-Nijenhuis structure. We shall assume that P and Q are such 
that at generic point of A^, the recursion tensor has m distinct (double) eigenvalues 
Aj. If moreover the m eigenvalues Aj are functionally independent when regarded as 
functions on A^, then it has been proved 0] that (at least locally) other m functions /i, 
exist such that: 

(i) the functions (Aj,/ij) form a system of coordinates on 

(ii) in this coordinate system, the Poisson tensor P is in canonical form, i.e. {Aj, /ij}p = 
5ij and {Aj, \j} p = {fii,fij}p = 0, and the recursion tensor is diagonal. 

The coordinates {Xi,fii) are called Darboux-Nijenhuis coordinates. The property (ii) 
completely determines the Q-Poisson brackets of the coordinates: {Aj,/ii}g = Aj and 
{Xi, fij}^ = {Xi, Xj}^ = {fii, fij}^ = for z 7^ j. 

Suppose hi to be a set of m hamiltonians, independent and in involution with respect 
to P and Q (although not necessarily generated by Magri-Lenard recursion). Falqui, 
Magri and Pedroni 00 have recently found an intrinsic coupling condition with the 
recursion operator, ensuring that all the functions hi are separable in the Darboux- 
Nijenhuis coordinates. 

In the case of a degenerate PQ structure, one can sometimes perform a reduction 
onto a symplectic leaf of P, by projection along appropriate transversal vectorfields. 
This allows one to compute explicitly not only the coordinates Aj (which are the roots 
of the characteristic polynomial of A^, or rather of its minimal polynomial), but also 
the other coordinates /Xj, as the values taken by a suitable polynomial p{X) after the 
substitutions A = Aj (for the construction of p{X), the exact statements and the proofs 
we refer the reader to ||24|| ). 

For a trihamiltonian structure {Po, Pi, P2) with a nondegenerate Poisson tensor Pq, 
one is naturally led to introduce two recursion operators, 

Ni = Pi- Po"' and N2 = P2- Po\ (3.1) 

so obtaining a compatible pair of Poisson-Nijenhuis structures. In this case, under appro- 
priate conditions it is possible to obtain trihamiltonian Darboux-Nijenhuis coordinates, 
having the most simple and natural property that one could imagine in this context: 

Proposition 3.3: Let Ni,N2 be the two tensors defined by ( |3.1|) from a compatible 
trihamiltonian structure on a 2m-dimensional manifold. If 

1. all the eigenspaces of Ni and N2 coincide (equivalently, Ni and N2 have the same 
centraliser) ; 

2. both Ni and N2 have m distinct eigenvalues, forming together a set of 2m inde- 
pendent functions, that we denote by Aj and fj,i respectively {i = 1, . . . ,m); 
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3. for any pair of eigenvalues Xi and Hi, corresponding to the same common eigenspace, 
one has 

then the eigenvalues of Ni and N2 (respectively denoted by Xi and Hi) form a Darboux- 
Nijenhuis coordinate system. 

Proof: For any Nijenhuis recursion operator and any of its eigenvaules A, it is 
always true (see [Q) that N*dX = XdX. The fact that A^^i has m independent eigenvalues 
implies that all its eigenspaces are bidimensional, and the same holds for N2. We have 
also assumed that the eigenvalues of Ni and N2 corresponding to the i-th common 
eigenspace are independent functions. Therefore, the eigenspace itself is spanned by the 
two differentials dXi and dfif 



N^dXi = XidXi, N^dXi = fiidX 

N^dni = Xid/j^i, N2djjLi = jJLidjjLi. 



From the definition ( |3.1| ) and the trivial fact that A'^iPq = Pi = PqN*, at any point 
where Xj 7^ one has 

{A,,A,}^ = ^{PodXi,XjdX,) = ^{PodXi,N*dXj) 

= ^{P,N*dXi,dX,) = ^{X„X,}^^ 

and, since Xi 7^ Xj for i 7^ j, one should have {Xi, Xj} ^ = for all In a similar way, 
using the recursion tensor N2 one obtains = for all Furthermore, 



1 1 

{Xi,fj.j} = —{PodXi,Xjdfij) = —{PodXi,Nldfj.j) 
Aj Aj 

= —{PoN*dXi,dij,j) = ^{Xi,iJ,j} 
Aj Aj 

which entails {Xi,fij}^ = for i 7^ j. The above results extend by continuity to points 
where Xj = 0. The additional normalisation condition (3) ensures that (Aj, are canon- 
ical coordinates for Pq; by construction, both tensors Ni and N2 are diagonal in these 
coordinates, which therefore are Darboux-Nijenhuis for the trihamiltonian structure. □ 

In other terms, if two Poisson pencils are available, and the two recursion op- 
erators are "independent" and "compatible" in the sense given above, then all the 
Darboux-Nijenhuis variables are obtained as eigenvalues (in the usual bihamiltonian 
setup, only half of these variables are defined as eigenvalues). The theorem also clarifies 
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that trihamiltonian structures generated by single recursion operator, Pi = NPq and 
P2 = NPi = N'^Pq, sometimes encountered in the hterature, are not suitable for our 
purposes. 

To obtain a twofold Poisson-Nijenhuis manifold also when the Poisson tensor Pq is 
degenerate, we exploit a projection technique, which has been already used for a different 
purpose in . 

Let us start with two compatible Poisson tensors Pq and Pi. The rank of Pq 
is assumed for simplicity to be constant on the phase manifold Ai (in the sequel, 
dim = 2m + k, k being the corank of Pq). The kernel of Pq is then pointwise spanned 
by the differentials of k Casimir functions c°. The main ingredient in our approach is a 
set of k independent vectorfields Za (spanning a distribution Z) having these properties: 

i) normalization: ^a(c^) = 5f ; 

a) integrability: [Za, Zp] G Z; (3.2) 

Hi) symmetry for Pq: Cz^Pq = 0. 

The normalization property (i) entails that all the vectorfields Z^ are transversal to the 
symplectic leaves of Pq. The concrete possibility of finding such a set of transversal 
vectorfields is not ensured a priori. Here we shall assume their existence, but in the last 
section of the article we will give an explicit recipe to find them for a relevant class of 
trihamiltonian structures. 

The vectorfields Za induce a decomposition on TM. and T*Ai: 



aT* KA / ^-L = {Za, 6')dc° 



MX e TM 



= X{c^)Za 

Xji = X — X_L 



Lemma 3.4: The decomposition satisfies the following properties: 

{d//,Za) = 0, 

x//{cn = 0. 

Proof: {9//, Za) = {e-e^, Za) = {6, Za) - {6, Zf,)Za{cf') = 0; 

X//(c°) = X(c") - Xx(c") = X(c°) - X(c'^)Z/3(c") = 0. □ 
Using the above-defined decomposition, one proves that 
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Lemma 3.5: The assumptions ( |3.4 ) imply [Za, Zp] = 0. 

Proof: [Za,Zp] E Z entails = 0, but one also has [Z^, = because 

\Z^,Z^\{c') = Z^{b})-Z^{bl) = ^. □ 

We can introduce a new tensor Pi by setting 

Pi(^^,0) = Pi(^^//,0//); (3.3) 

the kernel of the "deformed" tensor P\ contains all the differentials of the Casimir 
functions of Pq- Still, the tensor P\ is not automatically a Poisson tensor. 

From the beginning, we have assumed that C-z^P^ = (|3.2|) . This is equivalent to 
requiring that Pq be projectable along the flows of the vectorfields Z^, i.e. for any pair of 
functions (/,(?) such that Z^if) = Za{g) = everywhere, their Poisson bracket should 
also be constant along the same flows: Za{{f, g}p^^) = 0. 

So far, nothing ensures that the new tensor Pi is projectable in the same sense. 
Whenerver it turns out to be so, i.e. 

CzJi = 0, (3.4) 

then the Pi -bracket Pi{df,dg) can be reduced to Z^-invariant functions as well. For a 
Zcf-invariant function / one has df = {df) //, so the Pi -bracket of such functions coincide 
with their Pi-bracket. Then, if Pi is projectable, the Jacobi identity is straightfowardly 
proved for Za-invariant functions. We claim that (|3.4| ), supplemented by some auxiliary 
conditions, ensures both the Jacobi identity and the compatibility with Pq for the Pi- 
bracket of arbitrary functions on Ai. 

Proposition 3.6: If (i) all the vectorfields Z^ fulfill both and 

(a) all the Casimir functions c" are in mutual involution with respect to Pi, and 

(Hi) the functions generate bihamiltonian vectorfields, i.e. one can find k functions 

such that V" = Pidc" = Pod/i°, 
then Pi is a Poisson tensor compatible with Pq. 

The proof requires several steps. 

Lemma 3.7: If the vectorfields = Pidc" are hamiltonian also with respect to Pq, 
then the Poisson tensor Pi and the tensor Pi differ by a Lie derivative of Pq : 

Pi = Pi + Cx^^Pq, (3.5) 

where (summation is understood on repeated indices) 

Xp^ = h'^Za (3.6) 
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Proof: One has 



PM,dg) = P,{{df)i/,{dg),/) = Pr{df-{df)^,dg-{dg)^) 
= Pi(d/, dg) + Z^U)Zp{g)P^{dc^, dc") - 

-Z^{g)Zfsig)P^idf, dn - Z^{f)P^{dd', dg) 
= P,{df, dg) - Z^{g)Zp{g)Poidf, dh^) - Z^{f)P,{dh'^ , dg), 

which proves the statement since 

{Cx,^Po){df,dg) = h'^[Z^{{f,g}p^)-{ZM),9}p,-{f,Z^i9)}p,]- 

-ZM){h'',g}p,-Z^{g){f,h"}p^ 
= [{Cz^Po){df,dg)] - ZM){h",g}p, - ZMif^h'^yPo 

= -ZM){h'',g}p,-z^{g){f,h"}p^^. 

This also proves that, under the given hypotheses, Cxp^Po = V"' A Z^- □ 

If condition (iii) in the statement of Prop. (3.6) is fulfilled, then ( |3.4D can be rewritten 
in terms of the tensor Pi, using the previous Lemma: 



Lemma 3.8: Upon assuming ( |3.5| ), l \3.4 ) is equivalent to 



Cz^Pi = [Vl',Z^]AZp (3.7) 

Proof: Here and in the sequel we exploit the properties of the Schouten bracket, a 
bilinear operator on contravariant, antisymmetric tensors of arbitrary rank p (p- vectors). 
We refer the reader to [Vaisman] for the general theory; the properties that we shall use 
are the following: 

• if P is a p- vector and Q is a g- vector, then [P, Q]^ = (— 1)p^[(5, P]g; 

• if, moreover, P is a r-vector, [P, Q A R]^ = [P, Q]^aR+ (-l)«(P+i)g A [P, R]^; 

• if X is an ordinary vectorfield, then [X, P]g = CxP] 

• Cx[P, Q]s = l^xP, Q]s + [P, ^xQ]s; 

• if P is a bivector, then it is a Poisson tensor if and only if [P, P]g = 0; 

• if P and Q are both Poisson tensors, then they are compatible if and only if 

[p,g]g = o. 
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Proving the Lemma is now straightforward: 

= —^Zc^Xp^Pq 

= -[z^,v^]^Zp-v^ ^[z^,Zp] 

□ 

The condition ( ^.T] ) is also a projectabihty condition: it means that the tensor Pi 
itself can be reduced to Za-invariant functions, but in this case one cannot identify the 
image of Pi through projection (on the quotient space) with the restriction of Pi on a 
submanifold transversal to the vectorfields Za (for instance, a symplectic leaf of Pq). As 
a matter of fact, we have seen that the projection of Pi would instead coincide with the 
restriction of Pi. 

Now we can eventually prove the main statement. 

Proof: We need to show that [Pi,Pi]g = and [Po,Pi]g = 0. The second equality 
is easy to prove. By assumption, [Po,Po]g = [Pi,Pi]g = [Po,Pi]g = 0. Taking the Lie 
derivative of [Po,Po]g one finds that [Po,£x-Po]s = vectorfield X, and from 

(|3.5|) it follows that [Po,Pi]s = [PqiPi + ^Xp^Pq\^ = 0. To prove the first equality, ob- 
serve furthemore that the vectorfields V"" are assumed to be hamiltonian for Pi, hence 
[Pi, V% = Cv^Pi = Q. Then 

[Pi,Pi]s = [Pi,Pi]s + 2[Pi,l^"AZ„]g + [£xp,Po,^"AZjg 
= 2[Pi,y-]gAZ„-2l^-A[Pi,Zjg + 

[CxpPo, V\ A - A [CxpPo, ^a]g 

= 2F" A V"] AZp+ [y", A Zp]^ NZ^-V N [Z„ A Z^]g 
= 2\/" A [Z„ y^] A + [\/", yf"] A A + A Zp\ A - 

A [Z,, yf"] A - r° A A [Z,, Zf,\ 
= 2^ A V"] AZf.-V'^A [Z„ V^] AZ/s-V'A [Z^, V^] A Zp 
= 

□ 

We have seen that Pq and Pi are compatible and can both be reduced by projection 
along the flows of the vectorfields Z^ on any symplectic leaf of Pq. Then, on each 
leaf one can define a recursion tensor and look for Darboux-Nijenhuis coordinates; but 
there is no evidence that such coordinates can be extended to some neighborhood of the 
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symplectic leaf in Ai. However, a Nijenhuis tensor A^^i can be directly defined on the full 
manifold Ai in the following way. To any vectorfield X over Ai one can associate the 
vectorfield X// which is, by construction, everywhere tangent to the simplectic leaves of 
Pq. Consider now any one-form 9x such that X// = PqOx', we set 

A^iX = PiOx. (3.8) 

The tensor Ni is actually independent of the choice of 9x, because the latter is defined 
up to an element of the kernel of Pq, which is also the kernel of Pi. From the definition, it 
follows immediately that NiZ^ = 0, and that for any hamiltonian vectorfield Xj = P^df 
one has NiXj = Pidf . 

Proposition 3.9: Under the same hypotheses of the previous proposition, the tensor 
Ni defined by ( ^^ ) is a Nijenhuis tensor. 

Proof: We need to show that the Nijenhuis torsion tensor T/v of A^i vanishes: 

T^iX, Y) = [iViX, NiY] - Ni[NiX, Y] - Ni[X, N^Y] + N^[X, Y] = 

The Nijenhuis torsion tensor acts pointwise on vectors; thus, it is sufficient to show that 
T/v(X, y) vanishes on each pair of elements of a basis of the tangent space to A4 at 
any point. In our case, such a basis is provided by a set of Po^hamiltonian vectorfields 
spanning the tangent spaces to the symplectic leaves, and by the transversal vectorfields 
Za. We already know that NiZa = and [Za, Zp] = 0, thus TN{Za, Zp) = for any pair 
{a, (3). More generally, for any vectorfield X 

Tn{X,Z^) = Ni{Ni[X,Z^]-[NiX,Z^]) 
= NiiCzANiX) - NiCz^X) 
= NiiCzANi)Poex) 

but since Cz^Pi = j^z^Pq = 0, one has Cz^{Ni)Pq = 0. 

Finally, we evaluate the torsion tensor on two hamiltonian vectorfields Xj = P^df 
and Xg = Podg; for notational convenience, we set Yf = Pidf . This part of the proof is 
well known and can be found, in greater detail, in [Marsico]: it is shown there that the 
compatibility of Pi and Pq, i-e. [Po, Pi]s = 0, implies 

+ =X|^,,I^^ (3.9) 

Using the fact that Yf = NiXf, applying once again A^i to both sides of ( |3.9| ), and 
rearranging terms, one gets 

y{f,9U = NAXf,NiXg] + Ni[NiXj,Xg] - NnXf,Xg]; 
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on the oter hand, y{f,g}p = [^/,^] = Wi^ ^i^g] so 

T;v(X;,X,) = 0. 

□ 

These results suggest the following strategy: given a trihamiltonian structure, choose 
one of the Poisson tensors, say Pq] a complete set of Casimir functions c" can be directly 
read out from the "fundamental molecule", as well as the functions h°' and fc" used in 
(|3.6|) . In this case the vectorfields Pirfc" and P2dc'^ are automatically commuting, so 
the conditions (ii) and (in) of Prop. (3.6) are fulfilled. If one is able to find a complete 
set of vectorfields fulfilling all of (|3.2| ) and being symmetries of both Pi and P2 (in the 
sequel, we shortly write "a good set of transversal symmetries"), then one can apply the 
above procedure to each of the pairs (Pq; -Pi) and (Pq, P2) separately, obtaining in this 
way a pair of Nijenhuis tensors A^^i and N2 on the full manifold M.. 

This fact, however, does not yet ensure the existence of Darboux-Nijenhuis coor- 
dinates. The proof of the existence of Darboux-Nijenhuis coordinates has been given 
in Prop. (3.3) only for pairs of Nijenhuis tensors fulfilling some additional requirements, 
first of all that of being non-degenerate. We have not addressed the general problem 
of the existence of a canonical form for a non-regular Nijenhuis tensor; incidentally, to 
our knowledge non-regular Nijenhuis tensors have never been previously considered in 
connection with finite-dimensionl integrable systems. In the next section, we will rather 
characterise a particular class of trihamiltonian systems for which Darboux-Nijenhuis 
coordinates exist, are separation variables for the Hamiltonians occurring in the "funda- 
mental molecule", and can even be constructed without having to compute eigenvalues. 
Although the assumptions that we make could seem rather artificial and difficult to test 
in practice, we will eventually show that this class of systems is not empty, and contains 
relevant examples described by Lax equations (with spectral parameter). 

3.2 Sklyanin separation of trihamiltonian systems 

We first recall the separability criterion introduced by Sklyanin |T^. Let hi be m 
hamiltonian in involution for a nondegenerate Poisson tensor P, and let (Aj, /ij)i=i,...^m a 
system of canonical coordinates for P. If m functions Wi of m + 2 variables exist such 
that 

Wi{Xi,fii;hi{Xi,fii),...,hm{Xi,fii)) = 

; (3.10) 

W^m (Am,/im; /il(Ai,/ii), . . . , /lm(Ai,/ii)) = 

identically, then all hamiltonians hi are separable in the coordinates {Xi,fii). 
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This setting refers to the symplectic case (as it should be, since the notion of sep- 
arability of the Hamilton- Jacobi equation makes sense only in symplectic manifolds, 
namely on cotangent bundles). In the case of degenerate Poisson manifolds, one can dis- 
cuss separability on symplectic leaves. Upon reduction, one may expect to find Sklyanin 
functions Wi depending on auxiliary parameters labeling the symplectic leaves (in our 
case, the Casimir functions c"). This is indeed what happens in the example that we 
shall now discuss: first, we introduce a set of equations built using the common Casimir 
function /^^ (generating the "fundamental molecule") and the transversal vectorfields 
Za, and we list sufficient conditions for the roots of these equations to form, together 
with the Casimir functions c"), a system of generalised Darboux-Nijenhuis coordinates 
(Ai,/ii,c") on A4. Then, we produce m functions Wi, each depending on the i-th pair 
of coordinates {Xi,fii), on all the k Casimir functions c", and on other m arguments. 
Once the latter are replaced by the remaining m Hamiltonians h'^{Xj, fij, c") of the fun- 
damental molecule, the functions Wi vanish identically for any value of the coordinates 
c", thus on all symplectic leaves simultaneously. 

As previously, we assume that all the three Poison tensors (Pq^-PijA) have the 
same rank 2m, and the dimension of the manifold is 2m + k. We also assume that 
the common Casimir function fx^ is complete, i.e. that among its coefficients one can 
find m + k independent functions, including k Casimir functions for each of the three 
Poisson tensors (different tensors may indeed have some common Casimir functions). 
In the cases that we shall consider (for instance, the trihamiltonian spaces g[(r)"), the 
polynomial fxf^ has exactly m + k non-constant coefficients. 

The leading role in our construction is played by the derivatives of the common 
Casimir function /^^ along the transversal vectorfields Z^. We denote these k functions 
(still depending on the two parameters A,/x) by 

S^{X,fi) = ZMx,)- (3.11) 



The letter S is chosen because of the coincidence with Sklyanin's minors |T^, in the 
particular case discussed in section (4). 

Having assumed that the common Casimir function /a^ is polynomial in both pa- 
rameters, the functions S'„(A,/i) are polynomials as well. We shall prove that, whenever 
appropriate conditions are verified, the common roots (Aj, /ij) of the polynomials Sa{X, fi) 
are Darboux-Nijenhuis coordinates and fulfill Sklyanin's separability condition ( |3.11| ). 
Note that we use the abbreviated notation F{Xfi)\x.^^. for F{Xfi)\x=Xi,ii=Hi- 

Proposition 3.10: If a good set of transversal symmetries Z^ fulfills in addition the 
following requirements: 

1. all second directional derivatives of the complete common Casimir function fxfj, 
vanish identically, i.e. 

Z^iZf,{fx^))=0 for all a,/?; (3.12) 
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2. the polynomials S'q,(A,/x) = Za{f\fj) have 2m functionally independent common 
roots {Xi,fii}; 

3. the following equality holds, and for any i both sides are not identically vanishing 
for at least one pair {a, P): 



(3.13) 



then, the 2m functions (Aj, ^i) form a system of Darboux-Nijenhuis coordinates on each 
symplectic leaf of Pq, and moreover Za{Xi) = Za(fii) = 0. 

Proof: First, we compute the parameter-dependent vectorfield associated to the com- 
mon Casimir function fx^ under the deformed Poisson pencil (Pi — XPq): 



(A - APo) dfx, = (Pi + Cx.^Po - APo) dfx, = (Cx^^Po) df 



(3.14) 



But 



(Cxp^Po) dfx, = -J2 •^"(^' ■ Podh'' : (3.15) 

a=l 

in fact, from the definition of the deformation vectorfields (p.6|) , one sees that the vec- 
torfield {Cxp^ Pq) dfx, acts on an arbitrary function g as follows: 



^=0 



[{Cx.Po) dfx,] (g) = -Y, ZMx,){h\ 9},^ - E Z.(9){fx,. h^} 

a=l a=l 
k 

= -Y,S^{X,fi)-{Podh"){g). 

a=l 

Taking the derivative along Za of ( |3.14| ), one obtains 

k 

(Pi - APo) dS^iX, /i) = - 5^ Sp{X, /i) ■ Pod[Z^{h% 

(3=1 

(we used the fact that Z^ is a symmetry for both Pq and Pi). 

Since {Xi, Hi) is a pair of roots of Sa, one finds that [(Pi — APq) dSa{X,jji)] 
and that 



Q = d[Sa{Xi,lJLi)] = [dSa{X,fi)] 



+ 



'dX 
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dXi + 



then, 



dX 

On the other hand, 



{Pidfii - XiPodfii) = 



0, 



having assumed ^/^[^^(A, yu)] = and Sa{Xi, fii) = 0. Let S be the k x 2 matrix whose 



rows are 



' dSg dSg 



and let S(j) denote the same matrix after the substitution (A, fi) 



{Xi,fii). The condition ( p.l3| ) entails that S(i) has maximal rank for alH = 1, . . . , m so 
the previous results imply that 



Za{Xi) 

PidXi 
Pidiii 



XiPodXi 
XiPodfii, 



and the same holds for P2. Repeating now the same argument used in the proof of 
Prop. (3.3) one obtains that {Aj,Aj}^ = and = for all and that 

To get the remaining canonical bracket, one should instead rely on ( |3.13| ): for any 
pair of polynomials F(A,/i) and G(A,/i) such that F(Ai,/ij) = G(Aj,/ij) = 0, one has 



{F(A,/x),G(A,/i)}, 



fdFdG _dFdG\ 



then, if ( |3.13D holds and both sides are nonvanishing, one concludes 

{Ai,/Xi} = 1. 



□ 



Next, we prove Sklyanin's separability condition for the variables constructed ac- 
cording to the previous proposition. 

Proposition 3.11: Let (A,, /Xj, c")j=i,,,,,m,a=i,...,fc he the generalised Darboux-Nijenhuis 
coordinates associated to the projection along the vectorhelds according to Prop. (3.10) 
For each i = 1, . . . ,m 

hf^\x,,t,,=Pii>'hf^i) (3-16) 
where Pi{X, /i) are polynomials with constant coefEcients. Hence, the coordinates (Aj, /i,) 
and the remaining m hamiltonians (restricted to the symplectic leaf S ) fulfill the 
separability condition with 

Wi{Xi,fii;hl) = /amIa,,^., -Pi(Ai,/ii) 
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Proof: From the previous proposition, we know that 

^a(A„/ii)=0, Z,(A,) = 0, ZM = 0- (3.17) 

these imply that in the coordinate 



together with the normahzation condition 
system {Xi^Hi^c^) one has 



d 



(3.18) 



We denote the function obtained by replacing the spectral parameters (A, yu) with the 
pair of coordinates (Aj,/ii) by f(i) = fxfi\x.^.- On account of ( p.l7|) , 

thus the m functions /(j), which would in principle depend on all the coordinates 
{Xj, fij,c"), actually do not depend on the Casimir coordinates c". To prove that they 
depend only on the pair (Aj,/ij), we exploit the properties of Darboux-Nijenhuis coor- 
dinates: for any function g, one has 



Po 



and {gAj}. =^j^- 



Thus, if for some i one finds a function g such that {g, \j} _ — \i{g , \j} ^ =0 identically 

for any j, then (Aj — = and the function g cannot depend on fij for j ^ i. If 

furthermore {g,^j} — fii{g,fij} vanishes for any j as well, then g depends only on 

P2 ■''o 

the pair (Aj, fii). So what we need to prove is that, for any j ^ i, 



The differential of the function /(j) is given by 



0. 



(3.19) 



dfXi. 



Thus, the vectorfield defined by applying the tensor (Pi — AjPo) to the differential of 
the function /(«) can be obtained from ( p.l4| ) by replacing the parameters (A, /i) with 



the z-th pair of coordinates (Aj, /ij), and adding two terms proportional to and 
respectively. Applying this vectorfield to a coordinate Aj, with j ^ i, one gets: 



+ 



({/ii,A,}^^ - A,{/i„A,}^J ^ 



+ 
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By hypothesis, Xj is in involution with both Aj and fii for j ^ i\ therefore, only the 
first line survives, but due to (|3.14|) the r.h.s. is equal to — X]a=i '^'0(^4, Aj}^ , 

which once again vanishes on account of (|3.17|) : half of ( 3.19|) is proved. Repeating the 
whole argument for the coordinate /ij, upon replacing Pi with P2, one proves the full 
statement. □ 



Thus, the search for separation variables is completely translated into the problem 
of finding a set of "good" transversal vectorfields. We leave three questions open. First, 
suppose that one were able to find (in some other way) a set of generalised Darboux- 
Nijenhuis coordinates for the "deformed" Poisson triple (Pq, -Pi, -P2); would these coor- 
dinates be roots of the polynomials Sa (without imposing further conditions)? 

Second question, is requirement (3) in Prop. (3. 10) really necessary, or is it already 
implied by the previous assumptions? We could not find any concrete example in which 
( |3.2| ) and the requirements (1), (2) of Prop. (3. 10) are satisfied, but (3) fails to hold. 
One could then suspect that the equality (3) can be derived from the other (simpler) 
assumptions. The trouble with the requirement (3) is that it is both uneasy to check 
and lacking a clear geometrical significance^; nevertheless, we have not yet succeeded in 
replacing it with another condition equally ensuring that {Aj,/ij}^ = 1. 

Third open problem: under our assumptions, we have obtained Sklyanin's separation 
condition with Wi = f(^i) — pi{Xi, fii). In the algebro-geometric setting, one deals with a 
more particular situation, namely Pi{Xi, /ij) = p{Xi, fj,i) for a fixed polynomial p{X, fj,) not 
depending on i, so all separation variables are (pairwise) roots of a single polynomial 
fxfi ~p{^jf^)y defining the spectral curve of the system. At the moment, we do not 
know which additional conditions would ensure this stronger type of separability, which 
is encountered in the examples that we discuss below. In the next subsection we shall 
produce an example showing that the occurrence of a single spectral curve does not 
follow automatically from our assumptions. 



3.3 Canonical form of trihamiltonian structures in separation 
coordinates 

As we have seen, given a trihamiltonian structure (Po,Pi,P2) and a complete common 
Casimir polynomial fx^, one needs to find a good set of transversal vectorfields to pro- 
duce separation variables. Finding such vectorfields is, in general, a difficult task. On 
the other hand, upon assuming that such vectorfields exist, one can explicitly compute 

^In principle, it would be possible to consider the manifold x (with the spectral parameters 
A, ^ regarded as two additional real coordinates), endowed with the direct sum of the structure Pq 
on Ai with the Poisson structure on K^for which (A, /i) are canonical coordinates. Then, (3) could be 
rephrased by saying that all the Poisson brackets of the functions Sa (regarded as functions of 2m + k + 2 
variables) should vanish on each submanifold described by the equations X = Xi, fj, = fj,i. Still, this 
does not seem to help very much. 
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the components of all the relevant object as they would become in separation coordi- 
nates, as we shall see now; in some cases, this provides a concrete procedure to obtain 
"backwards" the change of variables. 

From the previous discussion, we know that in a generalised Darboux-Nijenhuis 
coordinate system: 

1. as far as the 2m x 2m block corresponding to the coordinates (Aj, /ij) is considered, 
the tensor Pq is in canonical form, while the tensors Pi and P2 are obtained from 
Po by applying the diagonal recursion operators Ni and N2, having the coordinates 
Xi and /ij, respectively, as (double) eigenvalues; 

2. the components of the complete tensors Pq, Pi and P2 in the coordinates (Aj, /ij, c") 
are obtained by simply adding k null rows and k null columns to the respective 
2m X 2m matrix (by "null" we mean that all the corresponding entries are van- 
ishing) ; 

3. the transversal vectorfields are coordinate vectorfields: Za = 

4. the relation between the original Poisson tensors Pi and P2 with the "deformed" 
ones. Pi and P2, is given by ( p.5|) ; 

5. each pair of conjugate coordinates (Aj, fii) is a root of the equation fx^—pi{\, /i) = 
for some polynomial pi. 

The latter information allows one to find the explicit expression of all the hamiltonians 
h^j as functions of (Ai,/ij,c"), once fixed the polynomials pi (which may be separately 
determined or arbitrarily chosen, as we explain below). In fact, let us assume as before 
that the common Casimir polynomial fx^ contains exactly m + k independent hamilto- 
nians h^j. Let us single out the k hamiltonians which are Casimir functions for Pq, which 
we denote as above by c", and denote the remaining (independent) hamiltonians as Ha, 
with A = 1, . . . , m. We impose the conditions 

= Pi(Ai,/ii) 

: (3.20) 

yUm) = Pm(A 

which form a linear system of m independent equations in the m unknowns Ha- Solving 
it, one finds Ha = hA{Xi, fJ^iyC"). 

Next, one produces the deformation vectorfields Xp-^ and Xp^ according to Lemma (3.7). 
Then, one can compute the components of the two Poisson tensors Pi and P2: 

Pi = Pi-/:xp,Po and P2 = P2 - Cx^^Pq. (3.21) 
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It is also important to remark that one can obtain as well the expression of the polyno- 
mials Sa{X,fJ.) = This fact will be used in the apphcations. 

To fix the ideas, we work out a concrete example. We remark that in this way 
we shall display a concrete case where the requirements of Prop. (3. 10) can be directly 
tested, showing that the set of systems fulfilling our requirements is indeed not empty 
(other examples can be produced in the same way, starting from different "fundamental 
molecules"). Take the 0t(3) "fundamental molecule" represented in ( |1.27| ). We set = 
/i2, = h\ and = h^; we also simplify the notation for the remaining hamiltonians 
by setting hi = Hq, h2 = h^, and = /ij. The common Casimir polynomial (apart 
from possible constant terms, which may anyhow be compensated in the polynomials 
Pi) becomes 

f\^l = hi + + ciA^ + /isyU + 02X^1 + csA/i^. (3.22) 

We leave the three polynomials pi, p2 and ps undetermined; for brevity, we write 
for pi(Ai,/ii), and so on. Imposing ( |3.20|) , the hamiltonians read as follows: 



hi = (A2/i3 - A3yU2)p(l) + (Aa/ii - \lli-i)p(2) + (A1/U2 - X2fJ'l)P{3) + 

[(A3 - A2)A2A3yUi + (Ai - A3)AiA3/i2 + (A2 - Ai)AiA2/i3] ci + 

[(/"a - 1^2)^2^2.1^1 + (/il - /U3)AlA3/i2 + (/i2 - /il)AiA2/i3] C2 + 
[(/W2 - /U3);U2/i3Ai + (/i3 - /il)/il/i3A2 + (/^i - 1^2) lJ'llJ'2>^z\ C3 
■ (Ai;U2 - A2/U1 + A2/U3 - A3/i2 + A3/ii - \lli-i)~^ 



(^2 - /i3)P(l) + (/is - ^^l)P{2) + - /i2)P(3) + 

[(A2' - A3')/il + (A3' - \^)^^2 + (Al' - \2^)^^,] Cl + 
[(A2 - Ai)/ii/i2 + (Ai - X-i)lilll^ + (A3 - A2)/i2/i3] C2 + 
[(/U2 - /Ui)/ii/i2 + (yUi - fi3)filfl3 + (/U3 - /i2)/i2/i3] C3 

■ (AiyU2 - A2yUl + A2Ai3 " A3Ai2 + hf^l - Al/i3)""^ 

(A3 - A2)P(1) + (A3 - A2)P(1) + (A3 - A2)P(1) + 

[(Ai - A2)AiA2 + (A3 - Ai)AiA3 + (A2 - A3)A2A3] ci + 

[(/il - /i2)AiA2 + (/i3 - /il)AlA3 + {fl2 - Ai3)A2A3] C2 + 
[(/Us' - /i2')Ai + (/il' - /i3')A2 + (/i2' - /il')A3] C3 

■ (AiyU2 - A2/il + A2/i3 - A3/i2 + A3/ii - Ai/i3)""^ 
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The deformation vectorfields ( |3.6[ ) are 



X 



Pi 



d_ 



_d_ 



Xp- 



P2 



d_ 



^277^ + ^3 



Starting from 



Pi 



( 



\ 
















M 













A3 









\ 
















-Ml 













-f 3 





















M3 









Cxp^Po and P2 



one can obtain the matrix expressions for the tensors Pi = Pi 

P2 — C-Xp^Pq- They coincide with the latter two matrices above, respectively, as far as the 
6x6 upper left blocks are concerned, while the remaining three rows and three columns 
are rather complicated for both tensors (and it would be pointless to write them down 
here). The fact that the tensors Pq, Pi and P2 above are pairwise compatible Poisson 
tensors, as well as the fact that the coordinate vecorfields = are symmetries of 
all of them, are trivially verified; then, is it enough to reverse the steps of the proof of 
Prop. (3.6) to show that Pi and P2 are both Poisson tensors compatible with Pq. One 



can then check directly that the function ( p.22| ) with the three hamiltonians computed 
above is a Casimir function for Pq, Pi and P2, as expected. 
The polynomials S'Q,(A,/i; Ai,/ij) are 



Si 



(A3 - A2)A2A3/ii + (Ai - A3)AiA3/i2 + (A2 - Ai)AiA2/i3 

Al/i2 ~ A2/il + A2/i3 — A3/i2 + As/ii — Ai/i3 
(A2' - A3')/Xi + (A3' - Ai2)^2 + (Al' - A2')/i3^ , , 



+ 



Ai/i2 - A2/il + A2/i3 - A3/i2 + A3yUi - Xi^l^ 

(Ai — A2)AiA2 + (A3 — Ai)AiA3 + (A2 — A3)A2A3 

AlAi2 - A2/il + A2Ai3 - A3Ai2 + A3yUi - AiyU3 



/i + A' 



^2 



(/^3 - A*2)A2A3/ii + (/ii - /i3)AiA3yU2 + {^^2 " ^^l)>^A2^^Z 



Al^2 — A2/il + A2/i3 ~ A3/i2 + A3/ii — Ai/i3 
(A2 - Ai)/ii/i2 + (Ai - A3)/ii/i3 + (A3 - A2)/i2/i3 

Ai/i2 — A2/il + A2/i3 — A3/i2 + A3^i — Ai/i3 
(/il - /X2)AiA2 + (/i3 - ;^l)AiA3 + (^^2 - ;^3)A2A3 

Ai/i2 — A2/il + A2/i3 — A3/i2 + A3^i — Ai/i3 



+ 



A + 



jj, + Xjj, 
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{1^2 - /•t3)/^2/^3Al + (/i3 - /^l)/^l/i3A2 + (/il - /i2)/il/^2A3 
A1/U2 - A2/U1 + A2/i3 - A3/i2 + As/ii - Ai/i3 
(//2 - /il)/il/^2 + (/^l - /^3)/^lA^3 + (/^3 - /^2)/^2/i3 

Ai/i2 - A2Aii + A2Ai3 - A3Ai2 + As/ii - Ai/is 

(;U3^ - /^2^)Ai + (/ii^ - /i3^)A2 + (/X2^ - /^1^)A3 
Ai/i2 — A2/il + A2/i3 — A3/i2 + As/ii — Ai/i3 



+ 



A + 



/i + /i 



This set of polynomials satisfies all the conditions of Prop. (3. 10): first, they do not 
depend any more on the coordinates c"; second, the reader can check straightforwardly 
that (Ai,/ii), (A2,/i2) and (A3, /is) are pairs of common roots of the polynomials Si, S2 
and 5*3; as far as the third condition is concerned, one should check that for any i = 1,2,3 



the two matrices 



{Sa, Si3}po 



and 



as„ asa 



ds„ as, 



(with a, /5 = 1, 2, 3) are 



ax a^j, afj, ax 

not identically vanishing and coincide. Direct computation shows that this is indeed the 
case: for i = 1, the two matrices are both equal to 







(A2-Ai)(A3-Ai) 



(Ai-A2)(A3-Ai) 

(Al-A2)(/U3-/ii) + 
\ +(/Ui-/i2)(A3-Ai) 



(A2-Al)(/i3-yUl)+ \ 

+ (/^2-/^i)(A3-Ai) 

(/^2-/il)(/^3-/^l) 



(/il-/i2)(/i3-An) 







/ 



and similar expressions, with the appropriate permutiations of indices, are found for 
z = 2,3. 

It is worthwhile to remark that one can produce, choosing arbitrarily the polynomi- 
als pi(A,/i), infinitely many families of hamiltonians which are separable according to 
Sklyanin's criterion, but do not coincide with the coefficients of a single spectral curve, 
unless one sets pi{\, /i) = P2(A, /i) = PsiX, jj) = p{X, /i). 

Notice that the choice of the constant polynomials Pi does not affect the polynomials 
5*0,. Actually, if one starts from a fixed trihamiltonian structure and is able to find 
vectorfields Za fulfilling all the requirement listed in Prop. (3. 10), then the polynomials 
Pi are determined a posteriori simply by plugging separately each pair of common roots 
of the polynomials Sa into fx^. The above reconstruction of the trihamiltonian structure 
goes in the reverse direction: the polynomials Pi are arbitrary and determine at the same 
time the hamiltonians and the Poisson tensors Pi and P2. 

Thus, in our framework the constant terms in Sklyanin's separation polynomials Wi, 
and a fortiori in the spectral curve (whenever it exists), are not directly encoded in the 
hamiltonian structure underlying a dynamical system and its symmetries (some aspects 
connected to the arbitrarity of the constant part of the the spectral curve equation have 
been addressed by J. Harnad in 123]). We have seen that the constant polynomials 
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Pi are determined by the choice of a set of transversal vectorfields, or equivalently of a 
system of separation variables: possible different sets of polynomials Pi - and, eventually, 
different spectral curves - are associated to different sets of separation variables. Indeed, 
from the example given above one might infer that different choices of the polynomials 
Pi lead to different hamiltonians (hence, to distinct dynamical systems), but in fact 
these are - by construction - nothing but the same hamiltonians expressed in two 
different coordinate systems: the same holds for the components of both Pi and P2. 
The transversal vectorfields, on the contrary, would have the same components in both 
coordinate systems, but the vectorfields themselves would not be the same, as happens 
for the deformed structures Pi and P2. 

In conclusion, the spectral curve appears to be an additional datum with respect to 
the purely hamiltonian structure of an integrable system; its hamiltonian interpretation, 
however, is deeply connected to the existence of particular canonical coordinates, as 
(separately) suggested by Sklyanin and Magri. For trihamiltonian structures (with a 
suitable projection onto a symplectic leaf), we have found a sound connection between 
separation coordinates and the vanishing of spectral polynomials. 

4 Separation coordinates for Lax equations with spec- 
tral parameter 

In this section we apply the techniques discussed so far to a particular class of dynamical 
systems, represented by Lax equations with spectral parameter. More precisely, we shall 
restrict to the following situation: 

1. the Lax operator is a r x r matrix polynomial of degree n in the spectral parameter, 
L(A) = AA" + MiA"~^ H h Mn, with a constant leading term A e g[(r); 

2. the constant matrix A should commute only with linear combinations of its powers 
(including A^ = 1). Equivalently, if A is diagonalisable, it should have distinct 
eigenvalues; more generally, the canonical Jordan form of A should not contain 
Jordan blocks proportional to the identity of dimension higher than one. This 
property is generic on 0[(r), but the requirement rules out some cases considered in 
the literature 0] 0|. A consequence of this requirement is that A may be nilpotent, 
with A^ = 0, but A^ should not vanish for k < r. 

3. the variable matrices Mi are generic matrices belonging to gl(r): we are not con- 
sidering restrictions to proper subalgebras, nor other types of reductions. 

This setting includes classical models such as the Euler-Poinsot top (for r = 3, n = 1) 
and the Lagrange top (for r = 3, n = 2), in the sense which has already been explained 
in the introduction: the classical models are properly embedded into larger systems. 
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but can be recovered by simple restriction (of a subset of the flows) on the appropriate 
invariant submanifold (namely, the subalgebra of antisymmetric matrices). Other impor- 
tant examples such as the Kovalewska top, or the Dubrovin-Novikov finite-dimensional 
reductions ||18[ of the Gel'fand-Dickey soliton hierarchies, are strictly related to the sys- 
tems that we are considering but cannot be directly obtained by simple restriction "a 
posteriori" . The application of our framework to the periodic Toda lattice or to other 
models without a constant leading term in the Lax matrix has not been investigated 
yet. 

There are two possible approaches to the construction of multihamiltonian structures 
for Lax equations of the type considered. Most authors regard them as dynamical 
systems on loop algebras Ql{r) = g[(r)((A)), and use the i?-matrix technique to define 
compatible Poisson brackets, which can be reduced to finite-dimensional quotient spaces 



identified with the linear spaces of fixed-order polynomials in A 

According to another approach, one considers the direct sum of n copies of the Lie 
algebra g[(r), defines a suitable Lie algebra structure on this vector space (different 
from the direct product structure), and an appropriate scalar product; in this way, 
one gets a natural Lie-Poisson bracket on Qi{r)^. The other Poisson structures are 
obtained by a deformation procedure, i.e. are defined as Lie derivatives of the Lie- 



Poisson tensor along suitable vectorfields |Tl|. In this approach, the dynamical variable 



is a n-ple of matrices (Mi, . . . , M„), while the fixed matrix A occurs in the definition of 
the deformation vectorfields which produce the Poisson structures; the Lax matrix L{X) 
itself arises as a by-product, in connection with the Hamilton equations with spectral 
parameter which naturally represent the trihamiltonian flows. The two approaches are 
substantially equivalent for our purposes. We shall follow here the second approach, but 
we stress that most of the definitions could be rephrased in the i?-matrix language. 

As anticipated in the introduction, we are not giving here all the proofs. The proof 
that the tensors Pq and Pi defined below are compatible Poisson tensors can be found. 



for instance, in |TT| or The only proof which is included concerns the fact that 

the characteristic determinant of the Lax matrix is the fundamental common Casimir 
polynomial for our trihamiltonian structure. 

4.1 AfRne Lie— Poisson pencils 

Let us consider the linear space 0t(r)" = 0" Ql{r)] we shall denote elements of this space 
by M = (Ml, . . . , Mn), with Mj G gl{r). The scalar product on gl{r) defined by ( |L5|) is 
extended "componentwise" to 0t(r)": 

n 

{A,B) = Y,Tr{A-B,). (4.1) 

1=1 

Using this scalar product, the gradient of a function / : Ql{r)"' — > M is again an element 
of Ql{r)"', that we denote by (Vi/, . . . , V^/). We define on 0((r)"' a first Lie-Poisson 
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structure P^"-* (depending on Mj in a strictly linear way) by setting 

n n 



(4.2) 



i=l 



k=i 



It is convenient to represent the Poisson tensor P*^"^ as a matrix of linear operators 
acting on the column vector (Vi/, . . . , Vn/): 



/ [Ml,-] [M2, ■] ■■■ [M„.i, ■] [M„, ■] \ 
[M2, ■] [Ma,-] ■■■ [Mn,-] 



p(n) 



[M„_i, ■ ] [M„, ■ ] 
V [Mn, ■] 













Following [|Tl[, from this first structure it is possible to obtain a sequence of other n affine 
Poisson structures P^"'~^\ . . . , P^^\ all mutually compatible, by the iterative formula 



-1) 



k + 1 



C^pC"^ k = n...l 



(4.3) 



where the components of the deformation vectors field X are affine functions, determined 
by a the fixed matrix A G gl(r): 



X 



n - 1 




V 



\ / Ml \ / \ 

+ 



0/ 



v 



V / 



The first pencil of the trihamiltonian structure that we shall consider is defined by the 
tensors P(i) and P(°) of the sequence; their expression is 



n— 1 i— 1 

{f^9}^,,, = (Vi(7,[V„_i/,A]) + 5^(v,(7,[V„_,/,A]+5^[V„W,Mfc 

i=2 k=l 

+(V„(7, [M„,V„/]) 



n i—1 

U'^^U) = (Vi^7,[V„/,A]) + 5^(v.(7,[V„„.+i/,A] + 5^[V„-,+,+i/,Mfc 

i=2 k=l 
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In matrix representation, 



p(i) 








A] 



.A] 
Ml] 



Mn-z] 
M„_2] 














(4.4) 



p(0) 












.A] 
MA 



The Poisson pencil P*-^^ — XP^'^^ is called the affine Lie-Poisson pencil on Ql{r). From 
the matrix representation given above, it is easy to see that any function /a = X] 
such that 



Vfc/A-AVfc+i/A = 
[AA- + MiA"-i + --- + M„,V„/a] = 



l,...,n 



1 



(4.5) 



is a Casimir function of the affine Lie-Poisson pencil. This observation leads to the 
definition of the Lax matrix L{X) = AA" + MiA"~^ + ■ ■ ■ + M„, and one immediately 
sees that the trace of any power of L{X) fulfills ( [4.51) . Each Casimir function of this type 
generate a bihamiltonian hierarchy according to the prescriptions ( |1.19| ) and ( |1.2CI| ). 
The corresponding Hamilton equations with spectral parameter are equivalent to Lax 
equations for L{X): the A;-th fiow of the hierarchy, using the same notation as in ( 1.19| ), 
is represented by 



L(A)= L(A),V„/, 



(fc) 



(4.6) 



In the sequel, we revert to the notation used in the previous part of the paper, 
setting Pq = P'^'^^ and Pi = P^^\ However, on gl{r)"' the third compatible structure 
P2, necessary to construct the appropriate second pencil, does not belong to the above- 
defined sequence of affine Lie-Poisson tensors: in particular, it has nothing to do with 
the structure P*^^^ , defined by (f4.3|) for n > 2 (for this reason we had to adopt a different 
notation). 

Our method to find the third Poisson tensor P2 still starts from the linear Lie-Poisson 
structure ( [4.2| ), but instead of applying a deformation vectorfield we exploit the Lax- 
Nijenhuis equation. This equation has been introduced in ||^ in connection with the 
following question: ''Given a Lax equation, and a Poisson structure P for which the 
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traces of the powers of the Lax matrix are in involution, does it exists a second compat- 
ible Poisson structure Q such that the same constants of motion are iteratively linked 
in a Magri-Lenard hierarchy?" It turns out that, if such a second Poisson structure 
exists, the two derivatives of the Lax matrix L along the two vectorfields Pdh and Qdh 
generated by any hamiltonian h are hnked by the following relation: 

CgdhL = ^Cpdh{L'^) + [L,a{dh)], 

for some matrix a algebraically depending on differential of the hamiltonian h. In some 
cases, this equation allows one to determine completely the second Poisson structure 
Q. This may happen if the Lax matrix is not generic, and in particular if L is assumed 
to have a fixed degree with respect to some grading in the Lie algebra; for instance, 
when L is tridiagonal (Toda), or is a polynomial of fixed degree in a spectral parameter 
(the case we are dealing with). Then, its square has usually a different degree, and 
the unknown element a occurring in the r.h.s. become determined by a compatibility 
requirement, i.e. its commutator with L should cancel exactly the terms of higher degree 
in the derivative of L^. We refer the reader to for a complete presentation of the 
method. 

We omit the details of the computation in our case: a key point is that one should 
plug in the Lax-Nijenhuis a slightly modified Lax matrix polynomial, namely the "con- 
voluted" polynomial L*{X) = A + Mi A + ■ ■ ■ + MnX"'- Then, starting from the Poisson 
tensor Pq, one finds the following new Poisson bracket: 

n n 

{f^9}p, = Y.i^^3.Y.MkVk-^+lfA-AVk-^+lfMk)+ (4.7) 

i=l k=i 
n n i—1 

Yl (^'3, Yl Yl MkVk+i-i+ifMi - MiVk+i-i+ifMk^ . 

i=2 k=i 1=1 

For n generic, writing down the representation of this Poisson tensor as a matrix of 
linear operators, analogous to ([4.4|) , would be rather cumbersome. The reader can easily 
figure out the general form from the representations of the Poisson tensors respectively 
corresponding to n = 1,2,3: 

for n = 1, Pa = M{-)A - A{-)M 



for n 



Mi{- )A- A{- )Mi M2{- )A- A{- )M2 
M2{-)A-A{-)M2 M2(- )Mi -Mi(- )M2 
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for n 



( 



\ 



A/i(. )A- A( - )Mi 

AfaC ■ )A- A( - )M2 
Aj-3(. )A- A{. )M3 



M2( ■ )A- A( ' )M2 

M3{ ■ )A- A{ - )M3 + 
M2( ■ )Mi - Mi( • )M2 

Af3(. )Afi - Afi{. )A.f3 



M3(- )A- A(- )A/3 
Af3( . )Mi - Afi( . )Af3 
Af3( ■ )M2 - M2(- )M3 



For n = 1 one recovers the Poisson structure of Morosi and Pizzocchero already 
described in the Introduction. The quadratic Poisson structures for n > 1 have never 
been presented in the previous hterature, to our knowledge. In the i?-matrix language, 
they can be obtained by a suitable modification of the so-called Sklyanin bracket pO^ 



4.2 Fundamental Casimir polynomial 

We have anticipated in the introduction that the characteristic determinant of the 
Lax matrix L{X) provides a complete common Casimir polynomial for the two pen- 
cils Pi — APq and P2 — /iPo- We shall now prove this statement. 

We already know that the traces of the powers of the Lax matrix are Casimir func- 
tions for the first pencil, as they fulfill ( |4.5| ). The same holds for the coefficients of each 
power of /i in the characteristic polynomial /a^^ = det(L(A) — /il), since these coefficients 
are functionally dependent on the traces of the powers of L{\). Then, what remains to 
check is that the polynomials in /i which occur as coefficients of each power of A in /a^ 
are Casimir functions for the second pencil. 

In the following computation, we denote the components of the gradient of 
by Vi = Vif\^^ (notice that each matrix still depends on both A and 11), and we 
systematically use the fact that 

Vk = XVk+u i.e. \4 = A^-V„. (4.8) 

The condition (P2 ~ f^Po) dfxfi = translates into the following system of n equations: 

' [fiV,, A] - ELiiMkVkA - AVuMk) = 

- YlU^kVu-^+iA - AVk-^+^Mk) 
. - ELi - MlVu+l-^-,lMk) = for 2 = 2, . . . , n 

which, with some algebraic manipulations taking account of ( |4.^ ) and of some of the 
equalities following from (Pi — APq) df\^ = 0, can be shown to be equivalent to 

(L(A) - VnMi = MiVn (L(A) - /il) for alH = 1 . . . , n (4.9) 

Now, for the characteristic polynomial fx^ = det(L(A) — /il) one has 

K = V„/a^ = /a^-(L(A)-/xI)-\ 
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therefore equation ([4 .91 ) is straigthforwardly satisfied. 

Under the hypotheses hsted at the beginning of this section (the phase space is 
the full space 0l(r)"', without constraints, and A is suitably generic), the non constant 
coefficients of the characteristic equation are functionally independent and define exactly 
|nr(r + 1) hamiltonians. Fig. 1 displays the corresponding "fundamental molecule". 
One sees directly from the diagram that, for each one of the Poisson tensors Pq, Pi 
and P2, the set of hamiltonians hj includes exactly nr Casimir functions, therefore the 
rank of the Poisson tensor cannot exceed nr^ — nr = nr{r — 1); on the other hand, 
the remaining ^nr{r — 1) hamiltonians are in mutual involution, so the rank of the 
Poisson tensor cannot be less than nr{r — 1). Therefore, the fundamental Casimir 
polynomial fx^ is complete. If A does not fulfill our requirements, what happens is that 
(i) the Poisson structures Pi and P2 have a larger kernel, so there are other independent 
Casimir functions not occurring in the characteristic determinant; (ii) some coefficients 
of the characteristic determinant vanish identically, and by consequence (iii) one cannot 
find properly normalised transversal vectorfields using the recipe presented below. 



4.3 Transversal vectorfields and separation coordinates 

The next step towards the construction of separation coordinates consists in finding a 
set of Po^transversal vectorfields, normalized on the nr Casimir functions which one gets 
from the fundamental Casimir polynomial. We state without proof the general recipe: 

1. Choose a matrix Wi G gi{r) of rank one such that: 

Tr{Wi) = 
Tr{WiA) = 

; (4.10) 

TriWiA"^-^) = 
TriWiA"^-^) = {-ly-^; 

this condition makes sense because we know that A'' for k < r (the r linear 
equations above do not determine uniquely Wi, but any such Wi will work). 

2. Compute the adjoint of the characteristic matrix L{X) — yul, that we shall denote 
by L"'"(A,/i); by definition, {L{X) — /xl) ■ L^{X,fi) = fx^l. The entries of L^{X,fj,) 
are the cofactors of L{X) — /iH, therefore V^{X,ix) is polynomial of order n(r — 1) 
in A and of order (r — 1) in /i. 

3. Let W\ be the matrix polynomial of the form 

Wx = WiX''-^ + W2X'^-'^ + --- + Wn. (4.11) 
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with 

W2 = iyi(Mo,2]I + Ml,2^ + --- + Mr-2,2^'"~^) 

; (4.12) 

The (r — 1) X (n — 1) coefficients scalar functions on gl(r)", which should 

be determined by the following condition: take the polynomial Tr(L"''(A, /i) VFa)/ 
for each power of /i separately, the n highest coefficients in A should vanish, except 
for the coefficient of A"^''~^-'+*^"^^\ which is always equal to one because of ( ^.lOp . 



Namely, the terms which should be canceled with the appropriate choice of Uij 
are the following: 

A*^, n(r-l) <k< 72(r - 1) + (n - 2), 

/^A^ (r2-l)(r-l) <k< (n - l)(r - 1) + - 1), 

/i"-lA^ <k< (n-l) 

This gives a linear system of equations for the coefficients Uij. The system is 
triangular and can always be solved. 

4. Denoting a tangent vector on 0[(r)" by v = [Mi, . . . , M„], introduce the n vector- 
fields 

v? = [O,---, 0,0,1^1] 
v° = [0,---,0,Wi,W2] 

\ (4.13) 

5. Take for /c = 1, . . . , r — 1 the product of Wx with the /c-th power of the Lax matrix 
L(A): 

H'f = (4.14) 

Consider the highest n terms in the expansion W^'^'^ = Wi'''^ X"''^ + '[¥2'^'^ X^^"^ + . . .+ 
Wn'^^X^ + . . . (in particular, one has Wi''^ = Wi A'^): other n{r — 1) vectorfields 
v*^, for z = 1 . . . , n and A; = 1, . . . , (r — 1), are defined analogously to (H:.13|) : 

= [0,---,0,0,<)] 
v^^ = [0,---,0,<),iyf] 

: (4.15) 
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For the nr vectorfields constructed in this way, the following holds: 

Proposition 4.12: The vectorfields are symmetries of Pq and fulfill the normaliza- 
tion condition 

v,'(/iU{r-fc-i)n) = S'^'Sjn^ for k, I = 0, . . . (r - 1) and j,m = 1, . . . ,n; (4.16) 
moreover, for all k, j, I and m in the given ranges, 

v,'K(W)=0 and [v^^v^]=0 (4.17) 

To ensure that the Po~ti'ansversal vectorfields provide a set of (separation) Darboux- 
Nijenhuis coordinates, one should also check conditions (2) and (3) of Prop. (3.7). How- 
ever, this verification may be postponed: let us tentatively assume that these conditions 
are satisfied. Then, one could find separation coordinates in which the Poisson tensors 
Pq, Pi and P2 assume the "canonical" form described in section (3.3); in particular, the 
affine Lie-Poisson tensor Pq becomes canonical in the usual sense. 

In separation coordinates, we already know the explicit form of the nr polynomials 
5*0 ( A, /i); the latter should coincide, up to the change of coordinates, with the polyno- 
mials Sj{X,fi) = \j ifxii) that one can also compute in terms of the original variables 
of the Lax equation. Hence, one should recover the separation coordinates by taking 
the common roots of the latter polynomials. Actually, it is sufficient to compute the 
common roots of the polynomials vj^ {f]iu)- The vectorfields vj^ are constant: namely, we 
have seen that they are defined by Mj = for i = l,...,(n — 1) and M„ = Wi A''. The 
corresponding polynomials Si{X,fi) are the derivatives of the characteristic determinant 
fxfj, along constant vectors fields, and therefore are nothing but linear combinations of 
cofactors of the characteristic matrix L{X) — /il, with constant coefficients determined 
by the matrices Wi and A. Therefore, we recover Sklyanin's recipe, with a definite 
prescription of the normalization to be used^ 

If one were able to compute ^nr{r — 1) pairs of common roots {Xi, fii) of the poly- 
nomials Si{X, fi), then one would only have to check a posteriori that the new variables 
{Xi, Hi) are canonical for the Poisson structure Pq. Since the polynomials S'f(A,/i) are 
independent linear combinations of cofactors of the characteristic matrix, {Xi, Hi) are 
also roots of the fundamental Casimir polynomial, and this would be enough to say 
that all the hamiltonians become separable in the new coordinates. In this case there is 
no ambiguity on the constant part of fxf^ itself, which turn out to depend only on the 
choice of A. Notice that in the new coordinates on 0l(r)"', given by the roots (Aj, fn) and 

^To be honest, exactly because Sklyanin's prescription does not fix the normalisation of the Baker- 
Akhiezer function, to prove the equivalence of the two procedures one cannot rely on the simple com- 
parison on the final results, but it would be necessary to compare each object involved in the two 
constructions, what we have not done yet. 
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by the nr Casimir functions of Pq? the vectorfields automatically become coordinate 
vectorfields, due to the normalization and to the property (|4.17|) . However, as we have 
stressed in the introduction, finding explicitly the roots of the polynomials 5'f(A,/i) is 
impossible in general. 

A more effective method is the following one: since we know the expression of the 
polynomials Sj{\,fi) in both coordinate systems, we can simply equate the correspond- 
ing coefficients for each polynomial to get a set of algebraic equations linking the two 
coordinate systems. Luckily enough, while the coefficients of Sj{X,fj,) are in general 
rather complicated rational functions of the separation coordinates, it is easy to see 
that a number of coefficients are linear functions of the entries of the Lax matrix L{X). 
It turns out that one can provide in this way a complete set of equations relating the two 
coordinate systems, which either are all linear (in the lower dimensional cases) or can 
be reduced to linear equations. In this way one can effectively compute the inverse co- 
ordinate transformation, i.e. express the original variables as functions of the separation 
coordinates. The mapping cannot be explicitly inverted in general, but it is sufficient 
for some purposes, for instance to check that the Poisson tensors transform as expected. 

Let us present a concrete example. We have already displayed in sect. (3.3) the 
form of the polynomials Sa for the system associated to the Lie algebra 0l(3). Let us 
now compute the same polynomials for the Lax matrix AX + M, with M G gl(3) and 
A diagonal with distinct eigenvalues (a,/3, 7) (this case corresponds to the generalised 
Euler-Poinsot rigid body, discussed in the introduction). We choose a set of orthonormal 
coordinates (xi, . . . ,0:9) in 01(3), setting 



M 



X-j 

-^{X4 — Xi] 



-j^{Xi + X4) 

Xs 

- X3) 



\ 75(^5 -X2) ^(X6 

The three Casimir functions for Pq (see ( p..27|) ) are 



hi 
hi 
hi 



/?7X7 + a'-yxs + a^Xg 
-{P + 'y)x7- {a + 'y)c 

Xt + Xs + Xg 



■j^{x2 + X5) ^ 
Xg J 



- (a + /3)xg 



(4.18) 



(4.19) 



Hence, there is a linear correspondence between the coordinates {xY,Xs,Xg) and the 
coordinates (ci, C2, C3). Choosing Wi to have all the rows equal to each other, following 
step (1) above one finds 



Wi 



1 




1 




1 


(a-/3)^(«-7) 


(7- 


-/3)^(q-/3) 


w- 


-7)^(0-7) 


(Q!-/3)(a-7) 
1 


(7- 


-/3)(a-/3) 
1 


(P- 


-7){a-7) 
1 



(4.20) 
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There are only three constant transversal vectorfields in this case, with M — Wi, M — 
WiA and M = WiA^ respectively. Taking the derivatives of fx^ = det{AX + M — 
//I) along these three vectorfields one easily computes the three polynomials Si{X,fi), 
S'2(A,^) and 6*3 (A, /i). Let us write down only the coefficients which are useful for our 
purpose: 



V2 



q; - 27 + (3)xi + {a-2(3 + 7)^2 + (2a - (3 - 7)^3 

+V2(7 - (3)x7 + V2{a - -f)xs + V2{P - a)xQ 

+ {a - (5)xi - {a- j)x5 + {(3 - j)x6 [(7 - (3){a - - 7)]~V 

+ ^ \j{a - 27 + P)xi - (3{a - 2/3 + 7)3:2 - a{2a - (3 - 7)2:3 
V2 L 

+^2(7 - /9)(-7 + a - /3)x7 - V2(a - 7)(a - /9 + 7)2:8 
+\/2(q! — /5)(q; — 7 + I3)xq + q;(7 — I3)xq — 7(0; — /9)x4 

- 7)^5] [(7 -/?)(«-/?)(«- 7)]"' A + .. . 

A// + ^= (2q;7 — /37 — q;/5)x2 + (q;7 — 2(5^ + Q;/9)a:3 

+ {-2aj3 + /97 + "7)a:i + V2(a - -i)(3xs - ^/2{a - f3)-fXg 
-a{-(3 + 7)0:6 + V2a{-f3 + 7)^7 + 7(0; - (3)x4 

-P(a--f)x5 [(7-/3)(a-/3)(a-7)]"V 
1 



H — ^ - (5{2a^ - (5^ - a(5)x2 
V2 L 

— q;(q;7 — 2/37 + '^l3)x^ — ^{—2af3 + f3^ + a'^)xi 
—\/2a^{a — 'y)xs + V2af3{a — f3)xg + a'^{—(3 + 7)2:6 
-a/27/?(-/? + 7)x7 - 7'(« - /9)a:4 
+/32(a - 7)0:5] [(7 - /5)(« - /?)(« - 7)]"' A + . . . 



V2 



{-f^a - (3^^ - a^(5 + a^7)a;2 + (7^" - h"^ + /3^a - /5^7)3:^3 

+(-/3^q; + /3^7 - q;^/3 + a^-i)xi - V2{-a(3 + a-f - (3-f)(a - 7)2:8 
—V2(—a(3 + 0:7 + /37)(q; - (3)xg — (a^ — + a(3){—(3 + ^)xe 
+V2{a^ - p-f + aP){-P + 7)2:7 + {-aP + a-f + p-f){a - P)x4 
+{—aP + 0:7 - P^) {a — 7)2:5 • 

•[(7-/3)(a-/3)(a-7)]-V 
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+ ^ [ - -i{-(3'^a + (3"^-^ - a^(3 + a^-f)xi - V2ajP{-P + 7)^7 

— /3(7^a — /57^ — + a^7)x2 — v^Q;/57(a — 'y)xs 
— 0(7^0; — /?7^ + — /5^7)x3 + V2al3'y{a — I3)xg 
— (—a/5 + 07 + Pj){a — l3)'~^Xi 

-{-a(3 + 07 - (3'^){a - 7)/3a;5 

+ (0:7 - /37 + a(3){-(3 + 7)q;X6 

•[(7-/3)(«-/?)(a-7r'A + ... 

Therefore, there are exactly six coefficients which are hnear in the coordinates Xj. One 
can check that they are independent and also independent from the other three equations 
already found from ( [4.19| ): the determinant of the full 9x9 linear system is identically 
equal to one. By equating the polynomials above to the corresponding expressions listed 
at the end of sect. (3.3), one finds the coordinate transformation. On the other hand, all 
three polynomials contain quadratic terms in A and /i, and to find the common roots of 
two of them one should solve an equation of order four in A (or in /x), which is possible 
in principle but gives a result which is of little practical use. 

The generalization of the Lagrange top can be treated in the same way. The phase 
space is gt(3)^, and as above we use as coordinates suitable orthonormal linear combi- 
nations of the entries of Mi and M2, that we denote by Xj, i = 1, . . . , 18. The constant 

(0 1 \ 
-1 . The Casimir functions for Pn oc- 
000/ 

curring in the characteristic determinant f\^j_ = det {AX^ + MiX + M2) are h^^, h^, h\, 
hi, Hq and hi; they are linear except for and which are quadratic. There are 6 
transversal vectorfields, three of which are constant. Each of the three corresponding 
polynomials, S2{X,fi), Sl{X, fi) and Sq^X,^!), has exactly three coefficients which are 
linear in the coordinates Xi, namely the coefficients of /i, A/z and A^. Then, one has 
a total of 13 independent linear equations which can be solved for the 13 coordinates 
Xi, . . . ,Xq, Xi5, . . . , Xis- The coefficients of A and A^ in the polynomials S^, S2 and 5*0, 
are linear in the remaining five coordinates Xio, . . . , X14, which can thus be computed as 
well. On the contrary, the direct procedure d la Sklyanin, i.e. computing the common 
roots of the bivariate polynomials Sa, is not viable due to the order of the polynomials 
themselves. 
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